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CHAPTER I 


General Introduction 


1.0 INTRODUCTION 

Modern civilization has become very complex with increasing 
demands of new and more resources in terms of materials, energy, 
habitats, etc. Rapid industrialization, fast population growth and 
associated pollution have created various kinds of problems 
leading to deterioration of our environment and ecology to a very 
undesirable level, lowering the carrying capacity of the earth to 
such an extent that if appropriate measures are not taken soon, 
the survival of the species will be threatened. Thus, conservation 
of resources, sustainable development and environmental protection 
by controlling the emission of pollutants are major 
socio-economic, scientific and technological challenges that our 
society faces today. 

In this thesis, therefore, the problems related to the fate 
of air pollutants and their removal in the atmosphere are studied 
using mathematical models. 

1.1 AIR POLLUTION 

The term environmental pollution may be defined as any 
undesirable addition to the environment which is harmful to man 
and his biosphere. In other words, it is the presence of 
undesirable contaminants in the environment in the form of gas, 
smoke, dust, fume, mist, waste disposal, etc. in such quantities 
and of such duration as to become injurious to human, plant and 
animal life, harmful to our living, to our enjoyment and to 
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historical monuments and other structures. Environmental pollution 
is generally divided into four broader classes : (i) Air Pollution 
(ii) Water Pollution (iu) Land/Soil Pollution, and Civ) Noise 
Pollution. The term social pollution is also sometimes used for 
abnormalities that arise in our society due to over population, 
crowding in habitats, rural migration, etc. Our main concern here 
is to study only certain aspects of air pollution using 
mathematical models. 

1.1.1 Air Pollutants 

Various gases and particulate matters emitted in the 
atmosphere by industrial stacks, automobiles, energy generating 
units, cooking appliances, etc., in concentrations significantly 
higher than their natural level, which are potentially harmful to 
man and his biosphere are referred to as Air Pollutants. The 
examples of air pollutants are particulate matters of both organic 
and inorganic compounds, smoke, oxides of carbon, sulfur and 
nitrogen and radioactive tracers, etc. Smoke is mainly a mixture 
of gases and unburnt carbon particles resulting from incomplete 
combustion. Significant amounts of air pollutants are produced by 
combustion of fuels and the type of pollutant depends on the 
nature of fuel used. For example, if coal is used, fly ash, sulfur 
dioxide and oxides of nitrogen are the major pollutants. It may be 
noted that sulfur dioxide is relatively short-lived in the 
atmosphere and gets readily converted to sulfur trioxide which is 
even more reactive and combines rapidly with any available 
moisture to form sulfuric acid mist which is the primary culprit 
in the phenomenon of Acid Rain. Similar phenomenon occurs with the 
oxides of nitrogen forming nitric acid as well. 
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Air pollutants may be categorized according to their nature 
at the time of emission from the source and the chemical reaction 
in the atmosphere, 

(i) Primary Pollutant : those directly emitted from a source or 
process, e.g., stack effluents, automobile exhaust. 

(ii) Secondary Pollutant : those formed in the atmosphere by means 
of chemical interactions among two or more primary pollutants or 
other atmospheric constituents, e.g., photochemical smog, 
formation of acid mist. 

Air pollutants are emitted in the atmosphere from various 
sources such as industrial stacks, power plants, automobile 
exhaust, burning of coal, fuel combustion, etc. These sources 
could be classified as point source (stacks), line source (highway 
vehicles) and area source (large industrial complexes). 

1.1.2 Effects of Air Pollution 

The effects of air pollution to human health, to ecology, to 
our property and to our enjoyment are well known. Air pollution 
can cause bronchial asthma, chronic bronchitis and other 
respiratory and heart diseases. Carcinogenic agents cause cancer, 
dust particles cause respiratory diseases and heavy metals like 
lead may cause poisoning. Carbon monoxide combines with the 
hemoglobin in the blood causing stress on those suffering from 
cardiovascular and pulmonary diseases. Radioactive fallout may 
cause biological effects like genetic damage or mutation. The 
health effects on human beings, include eye irritation, nose and 
throat irritation and increase in the mortality and morbidity 
rate. Some specific pollutants and their effects on human health 


are 1 isted below , 
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(i) Sulfur dioxide - an irritant, gas when inhaled affects the 
mucous membrane. Under certain conditions it oxidizes to sulfur 
trioxide which is much stronger irritant. 

Cii) Carbon monoxide - has strong affinity for combining with 
the hemoglobin molecules of the blood forming carboxyhemog 1 obin 
(COHb) which reduces the overall capacity of the blood to carry 
oxygen to the cells. 

(iii) Lead - emitted mainly from automobiles in urban 

atmosphere. Its effects include gastro-intest inal damage, liver 

and kidney damage, abnormalities in fertility and pregnancy and 
mental development of children. 

(iv) Insecticides - can affect the central nervous system and 

may attack other vital organs. Undesirable effects of indoor 

spraying may affect domestic livestock. 

(v) Hydrocarbon vapors - formaldehyde is primarily irritating 

and contribute to eye and respiratory irritation caused by 

photochemical smog. 

(vi) Radioactive isotopes - the major sources of radioactive 
air pollutants are nuclear reactors, testing of nuclear bombs in 
the atmosphere. The serious health effects are anemia, leukemia, 
and cancer, genetic defects and sterility as well as embryo 
defects and congenital malformations. 

Human health and ecology are also affected by toxic gases 
released in the atmosphere due to accident and one of the most 
disastrous episode in the industrial history was that of Methyl 
isocynate (MIC) gas leakage in a pesticide plant in Bhopal, 1984, 
India in which large number of people were reported to have 
suffered in different forms. Methyl isocynate decomposes easily 
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with moisture forming methylamine, and is highly toxic and 
unstable. (Fieser and Fieser, 1961) 

CH„NCO + H„0 > CH„N = C > CH„NHCO„H > CH„NH„ + CO„ 

o .c; J ^OH ^ 

MIC Methylamine 

MIC affects all living beings including vegetation. In a report 
(Arora colony, 1984), the following types of effects were noticed 
when some people were exposed for 1-5 min., 

0.4 ppm no effects 

2.0 ppm lachrymat ion , irritation in nose and throat 

4.0 ppm symptoms of irritation more marked 

21.0 ppm unbearable irritation of nose, eyes and throat 

Prolonged exposure to MIC destroys proteins and lipids in the 
lungs which leads to changes in the permeability of the lung 
membrane, getting filled with water causing death. The maximum 
permissible limit of exposure for industrial workers is only 0.2 
ppm for MIC, (Bretheric, 1981) 

Nuclear radiations from various sources constantly pose an 
environmental risk. A deadly mishap occurred due to nuclear device 
failure in Chernobyl, 1986, Russia, emitting huge amounts of 
radioactive material. This incident devastated a vast area around 
Chernobyl and gave rise to a nuclear cloud which remained cast 
over a large part of Europe for several days. In the case of a 
nuclear accident, dangerous radiations are emitted along with 
other harmful materials remaining in the atmosphere for a long 
time posing a threat to mankind. For example, an isotope of 
Strontium CStrontium-90) causes blood cancer and decays very 
slowly, Cesium-137 binds itself to soil and causes genetic changes 
and Iodine-131 has a tendency to cause thyroid cancer. 



Air 


P^Hution has also 


vegetation int 

Smog, dust, et 


adverse effect on plants and 
erfering with the plant growth and photosynthesis. 


C-. reduce the amount of light reaching the plants 
and cloggij^g tK^ 

^ stomati 


leaves of the 
the form of 


a reducing the carbon dioxide intake. The 
plants are severely affected by the pollutants in 
ci'osis, chlorosis, etc. The acid rain which is a 

major Pollutant 

affects forests, vegetation, crops, monuments and 

C8.LIS0S ^ C j. ^ 

ioation of water resources such as lakes, etc. 
destroying , 

P species and other marine ecosystems. 

Aip 

ion is also associated with the damage of property 

as well as oth 

economic losses. For example, due to pollutants, 

metal ^^r'rodes 

I'ubber cracks and looses its elasticity, leather 

becomes ^^ittl 

Sls-ss is etched, dyed materials may fade, etc. 

The damage + 

materials are done mainly through abrasion, 

deposition on 

ni'faces and by direct or indirect chemical attack. 

The Taj Mahal 

a miracle in marble, is also facing serious threat 

due to Commissi 

Pning of Mathura Oil Refinery about 40 kms. from 

the monument. 0^.,^ 

ifur dioxide mixes with atmospheric moisture and 
gets pppverted + 

o sulfuric acid which reacts with the marble 
(calcium ca^bonat«^ 

) resulting m corrosion and decoloration of the 

monument . 

VisiW*, • 

^'ISlbllity 


causing accident 
economic losses. 


may also be reduced due to formation of smog 
s, delay of flights, trains, etc. leading to 


1.2 dispersion OF 


AIR POLLUTANTS 


The 


di 


®Persion of air 


pollutants in the atmosphere is 


governed by th 

® processes of molecular diffusion and convection 

and is Riainly 

fected by various meteorological conditions such 



7 


as wind, temperature inversion, foggy atmosphere, topography of 
the terrain, etc. and types and number of sources. This process is 
also affected by several removal mechanisms such as chemical 
transformation, dry deposition, greenbelt plantation and 
precipitation scavenging, etc. In the study of pollutant 
dispersion, we are generally interested to find the concentration 
of the pollutant, after being emitted from its source, at a 
specified location and time. The concentration depends upon number 
of sources, their locations, emission rates and other 
meteorological and atmospheric conditions. Pasqu i 1 1 ( 1 962 ) , 
Ledbetter(1972) , DobbinsC 1979) , Wark and Warner ( 1 98 1 ) , 
SeinfeldC 1986) , Shukla et a/. (1987). 

The hot plume from the stacks goes up to a certain height, 
called plume rise, into the atmosphere, due to buoyancy and 
momentum and then spreads in both horizontal and vertical 
directions by the dispersion process. This plume rise depends upon 
many factors like exit velocity, temperature of the plume, wind 
velocity and diameter of the stack, etc. The height of the stack 
is also an important factor in determining the level of pollution 
at a given location. The total sum of the actual stack height and 
the plume rise beyond the stack exit is termed as effective stack 
height . 

1.2.1 Diffusion Equation 

Diffusion is the process by which matter is transported from 
one location in the atmosphere to another as a result of random 
molecular motions. The mathematical theory of diffusion in an 
isotropic media is based on Pick’s law of diffusion which states 
that the rate of transfer of diffusing material through unit area 
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of cross section is proportional to the concentration gradient 
measured normal to the section, (Crank, 1975) 

F = - D (1.1) 

ox 

where F is the rate of transfer of diffusing material per unit 
area of cross section, C is the concentration of diffusing 
material and x is the space coordinate. Here D, the constant of 
proportionality, called the diffusion coefficient, depends upon 
many factors including the concentration. The diffusion equation 
in isotropic media may be written as, (in vectorial notation) 

^ = div(D grade) (1.2) 

where D may be function of x,y, 2 , the space coordinates, and C. 

1.2.2 Convective Diffusion Equation 

Assuming continuity in the atmosphere and noting the species 
balance, the convective diffusion equation or atmospheric 
transport equation can be derived as. 



i = 1,2,3. . . (1.3) 


where x,y,z are cartesian co-ordinates, t is the time, u,v,w are 

the components of the wind velocity in x, y and z directions 

respectively and K , K and K are turbulent diffusion 

X y z 

coefficients. The terms and represent source (sink) and 

reaction (removal) terms, Dobbins( 1979) , Seinfeld(1986) . The 
diffusion coefficients are, in general, functions of space 
variables and concentration, however, they may be taken as 
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constants for simplification. The source term can also be included 
in the boundary conditions rather than in the equation itself. 

The above equation is widely used to study the dispersion of 
air pollutants with number of approximations and simplifications. 

1 . 3 REMOVAL MECHANISMS 

Various removal mechanisms are known to occur in the 
atmosphere which are effective in removing pollutants. Since the 
dispersal process is significantly affected by the removal 
mechanisms, it is absolutely essential to study the effects of 
removal mechanisms on the dispersion of air pollutants in the 
atmosphere. Some of the important removal processes are described 
below, 

Ci) Chemical transformation (reaction) - Many chemical reactions 
are known to occur in the atmosphere. Some of these may be very 
complex as that in the formation of photochemical smog. Others may 
be less complex as transformation of sulfur dioxide to acid rain. 
These chemical reactions are generally considered to be of first 
order and reduce the concentration of primary pollutants in the 
atmosphere. This mechanism can be used in removing the toxic 
species present in the atmosphere by introducing suitable liquid 
or gases which could react chemically with toxicants and 
neutralize their effects. 

(ii) Dry deposition - The term dry deposition means direct 
collection of gaseous and particulate matter on earth’s surface, 
soil or vegetative canopies. The deposition of pollutant species 
plays an important role in limiting the atmospheric concentration 
of the species. The vertical downward flux of a species is 
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represented in terms of an empirical parameter, called the 
deposition velocity , multiplied by the concentration of the 
material at some height i.e., 

Flux = v^C(x,y ,Zj , t) Cl. 4) 

where C is the concentration of the pollutant in the atmosphere. 
Modifying equation Cl.l) for turbulent diffusion coefficient and 
using equation (1.4) one can write, 

li = v^C(x,y, 2 ,t) at 2 = z, (1.5) 

This formalism applies to the diffusion of gases and fine 
particles when the gravitational settling velocity is negligible. 
The deposition velocity is certainly not a constant parameter and 
is subject to variability depending on the particular species or 
material that is being removed, on the meteorological parameters 
characterizing the state of surface layer, and on the nature of 
surface itself. Here v^ is a parameter that is proportional to the 
degree of absorptivity of the surface. At low emission heights, 
the removal of the pollutant by dry deposition must be taken into 
account in modelling process because under stable atmospheric 
conditions, the plume is restricted to a narrow boundary layer and 
the interaction of surface with plume can cause considerable mass 
removal . 

The interaction between the diffusing material and the 
surface is mainly due to total absorption (i.e., if earth is a 
perfect absorber, the concentration of material at z = 0 is zero) 
total reflection, and partial absorption where some particles are 
reflected and some are absorbed. The case of partial absorption at 
z = 0 is given by, (as in equation 1.5) 



at z 


0 


K = V C 

Z dz d 


( 1 . 6 ) 


and for total reflection = 0 and = a for total absorption. 

If the effect of gravitational settling of pollutants 
(particulate matter), which, however, can be neglected for gases 
and particles smaller than approximately 10 /im radius compared to 
turbulent diffusion, (Horst, 1984), is also to be included, the 
equation (1.5) modifies to 


K ^ + w C = V ,C at z = 0 

Z oZ S Q 


(1.7) 


where w is the settling velocity of pollutants, 
s 

If there is an impermeable layer at z = H which inhibits 
vertical mixing, called the inversion layer, the boundary 
condition at z = H can be deduced from equation (1.5) as. 


K |£ = 0 
z dz 


at z = H 


( 1 . 8 ) 


that is, when the effect of inversion is to be included, the 
condition at z — > oo is simply replaced by a condition of zero flux 
in the z-direction at the inversion height provided there is no 
leakage in the inversion base. 


(iii) Removal by greenbelt - The pollutants in the atmosphere 
can also be removed by suitable greenbelt plantation as plants 
have the capability of absorbing the pollutants. It may be 
mentioned here that the deposition, absorption characteristics of 
the greenbelt depend not only upon the distribution, size, height, 
nature of the canopy, etc. but also on the type of the pollutant, 
species of trees, shrubs, herbs, etc. planted in the greenbelt. 
Thus, the choice of the plantation should be such that the species 
planted be relatively resistant to the exposure of particular 
pollutant. The greenbelt acts like a sink for airborne material 
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because of its large surface area, the foliage of vegetative 
canopy serving as an effective receptor. If a greenbelt is planted 
suitably between the source and the region (i.e., a habitat, a 
historical monument, etc.) to be protected, it is possible to 
minimize the effect of pollution, the height of the source, 

thickness and the mean height of the greenbelt, determining the 
optimal location of the greenbelt. 

(iv) Wet removal or precipitation scavenging - Precipitation 
scavenging (the process in which atmospheric gases/particulate 
matters are absorhed/trapped in rain droplets falling on the 
ground) provides important removal mechanism of pollutants from 
the atmosphere. This mechanism is sometimes subdivided into 
rainout (i.e., the process within clouds where plume expands 

freely into the cloud layer and material is ’removed’ by getting 
dissolved in water or if it consists of particles, by getting 
entrapped in water droplets causing lowering of pollutant 

concentration) and washout (i.e., the process taking place below 
the cloud base for removal of airborne material by rain). 
Precipitation scavenging, not only acts like an important sink for 
trace gases in the atmospher but also influences the local, 
regional and global distribution of air pollutants in the 

atmosphere. 

1.4 PATCHINESS IN THE ATMOSPHERE 

It is pointed out here that the dispersal and removal 
processes are also affected by the ecological conditions on the 
ground and meteorological conditions in the atmosphere. For 
example, presence of a forest in a region would have different 
deposition rate affecting the dispersal process in comparison to 
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plains having little or no vegetation. Similarly these processes 
will be affected if rain is taking place in a region or it is 
covered by fog. Thus in real situations, the atmosphere may have 
patchiness with different dispersal, deposition and removal 
parameters in different patches and the modelling process must 
take this aspect into account. 

1 . 5 OUTLINE OF THE RELEVANT LITERATURE : MODELLING THE FATE OF 
AIR POLLUTANTS IN THE ATMOSPHERE 

As pointed out earlier, the dispersion of air pollutants is 
affected by meteorological conditions in the atmosphere, nature of 
terrain downwind from the stack, physical and chemical properties 
of the pollutant, number and location of stacks, etc. Several 
analytical studies have been conducted with some of the above 
mentioned factors but none of them accounts for all. 


1.5.1 Simple models 

The most comprehensive approach to transport theory is based 
on eddy diffusion model given by equation (1.3). The equation is 
quite complicated to solve analytically. It can be -simpl if ied to 
the following form. 


at “ax “ay "ax N ^^2 ^ 


(1.9) 


by assuming constant diffusivity coefficients, no chemical reactioi 
and other removal mechanisms in the atmosphere. 

Using following assumptions, 

(i) The process is steady state, i.e., dC/dt = 0 

(ii) The wind speed u is constant at any point x,y,z 

(iii) Major transport direction due to wind is along the x axis 
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(iv) The downwind diffusion is negligible in comparison to 
advect i on 

equation (1.9) can be further simplified as, 


,, ac _ a^c . ^ a^c 
ax ^,2 ^2 ^_2 ' 


K / K 
y z 


( 1 . 10 ) 


ay az' 

The approximate solution to equation (1.10) under the following 
boundary conditions has been given by Lowry and Boube 1 ( 1 967 ) , 

(i) C — > 00 as X — » 0 i.e., concentration is large at the poin 

source 


(ii) C — > 0 as x,y,z — > 
large distances from the 

(iii) ♦ 0 as z— ♦ 0 
z az 

(iv) 

0 


I I 

-d> 


uC(x,y,z) dydz 


00 i.e., concentration approaches zero 

source 

i.e., no deposition on the surface z = 

= Q ,x >0 i.e., rate of transport of 


at 


0 


pollutant downwind is constant and equal to the emission rate Q of 
the pollutant at the source, as follows: 


C(x , y , z ) = 


Q 


4irr(K K ) 
y z 


TT2 


exp 


r ^ 1 


( 1 . 11 ) 


, 2 ^ 2 ^ 2 . 1/2 

r=(x +y +z) 


The center line concentration at the ground level is given by, 


C(x,0,0) = 2 (1.12) 

4nx(KyK^)^^^ 

which is inversely proportional to x and independent of wind speed 

u, the two departures when compared with the experimental evidence 

176 

which shows that C is inversely proportional to (ux ), (Wark 
and Warner, 1981). Equation (1.11) shows that away from the centeri 
line the concentration decays exponentially in both y and z 
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directions. This mathematically means that C in the cross wind an( 
vertical directions may be "normally" distributed. Also th« 
decrease in the value of C in the x-direction is largely dependent 

on K and K , showing strong dependence on diffusion coefficients. 

y 2 . 

The other most widely accepted model which depends on the 
diffusion coefficients and exhibit the "normal" distribution as ir 
equation (1.11) is the Gaussian plume model. The Gaussian plume 
model solution for a point source at the ground level in absence 
of any removal mechanisms is given by, 


C( X , y , z ) = 


Q 


TtUOyO^ 


exp 


where Q is the source 


and o , a are the 
y z 

distributions related 



strength i.e. , mass 
standard deviations 
to mass diffusivities 


Cl. 13) 

emitted per unit time 
of the concentration 
K and K as , 

y z 


a 


2 

y 


2KyX 


2K X 
z 

u 


The equation (1.13) is further modified to, if emission of a 

gaseous pollutant takes place from an elevated point source at an 

effective height h , without reflection at the ground 

s 


C(x,y,z) = 


Q 

2itno^o^ 


exp 


[ - ^(4 ^' 4 ’ ]] 


(1.14) 


However, for a reflection at the ground, it takes the| 
following form, 

O A _ 2 r f -(Z-h )^^ r -(Z + h )^>-] 

Hi -^7^ J * ^7^ JJ 


(1.15) 


The above Gaussian plume equations, for non-reacting and 
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non-depositing pollutants from a continuous point source, in spite 
of their limited validity, have been extensively used for aii 
pollution studies. 

1.5.2 General models 

Various investigations have been made to study the dispersior 
of air pollutants in the atmosphere under varying meteorological 
conditions such as wind, temperature inversion using both eddy 
diffusion models and Gaussian models. Smith(1957), Turner ( 1 970 ) , 
Peters and K1 inzingC 1971 ) , Munn and Bolin(1971), Heines and 
PetersC 1973a) , Lamb and Se inf e 1 d( 1 973 ) , Nunge ( 1 974a , b ) , Hameed 
(1974a, b), Yeh and Huang(1975), Ragland and Denni s ( 1 975 ) , Sander 
and Seinfeld(1976) , Mehta and Ba 1 asubramanyamC 1 977 , 1 978 ) , Demultt 
(1978), Kao and Yeh(1979), Lupini and Tirabass i ( 1979 , 1983) , 
Draxler( 1980) , Nieuwstadt ( 1980) , Nieuwstadt and De Haan(1981), 
Wi 1 son( 1981 ) , LI ewelyn( 1982) , Karamchandani and Peters( 1983) , 
Robson( 1983) , B1 ondin( 1984) , Peterson( 1985) , Lor imer ( 1 986 ) , 
Lorimer and Ross(1986), Okamoto and Shiozawa ( 1 987 ) , 
Overcamp( 1988) , Lehmann( 1989) , Mihail et a7.(1989), Kashibhatla 
and Peters( 1990) , Singh et ai .( 1990a , b) , Eppel et a7.(1991), 
Runca(1992), Pandey et ai.(1992). In particular, Heines and Peters 
(1973a) presented an analytical model for continuously emitting 
infinite line and point sources with temperature inversion and 
constant wind velocity. 

In general, the wind velocity and dif fusivities are functions 
of space variables and time. When horizontal homogeneity is 
assumed, the wind velocity and diffusion coefficients may be taken 
as functions of height only. For general functional forms of these 
variables, no exact solutions of three dimensional diffusion 
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equation with chemical kinetics and other removal mechanisms are 
known. In a few cases when wind velocity and d i f f us i vi t i es are 
approximated, some solutions are available for two dimensional 
steady state equations. Smith(1957), Lebedeff and Hameed(1975) and 
Mehta and BalasubramanyamC 1978) have studied the problem of steady 
state two dimensional advective diffusive transport of a 
non-reactive pollutant emitted from ground based area sources with 
variable wind and turbulent diffusivity approximated by power law 
profiles with no flux condition at the boundary of the 
contaminated layer. Lupini and Tirabassi (1979) have proposed a 
’modified’ Gaussian model for two dimensional steady state 
advective diffusion equation for a non-reactive pollutant emitted 
from an elevated point source and also studied the effect of plume 
trapping beneath an inversion by this ’modified’ model. 
NieuwstadtC 1980) gave an analytical solution of the time dependent 
one-dimensional diffusion equation with more realistic diffusivity 
profiles than the power law profiles under stationary conditions 
for a non-reactive pollutant from an instantaneous area source 
assuming horizontal homogeneity in the atmosphere. Nieuwstadt and 
De Haan(1981) extended the above solution to non-stationary 
conditions in a developing boundary layer and applied the study to 
the special case of inversion-rise fumigation. 

It may be noted here that the Gaussian puff models are useful 
to simulate the instantaneous releases whereas Gaussian plume 
models are used for continuous releases whose duration is 
significantly longer than the travel time to a receptor. 
WilsonCl981) has given extended version of the Gaussian plume 
model to account for ensemble-averaged along-wind diffusion of 
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time dependent source release rates and presented explicit 
solutions for a suddenly started release with exponential decay 
rate. Using moments method, Lupini and Ti rabass i ( 1 983 ) obtained 
approximate solutions of the advective diffusion equation for 
horizontally homogeneous atmospheric boundary layer for a passive 
pollutant emitted from an elevated instantaneous point source. 
Overcamp (1988) developed models to simulate transient releases 
with exponential decay, of a pollutant that undergoes first order 
decay. The models developed are based on one-dimensional diffusion 
in the downwind direction coupled with the conventional Gaussian 
model in the vertical and cross wind direction. 

1.5.3 Effects of chemical reaction 

As we have stated earlier in section 1.3 that many removal 
processes are known to exist in the atmosphere, the atmospheric 
dispersion models should, therefore, be modified to account for 
these removal mechanisms for more realistic dispersion phenomena. 
A few mathematical models are available that takes into account 
the effect of chemical transformation on dispersion of air 
pollutants e.g. , Heines and Peters(1973b) , Nunge and 
VaidyanathanC 1977) , Peters and RichardsC 1977) , KewelyC 1978) , 
Naresh and Nath(1990). Heines and Peters ( 1973b) have analytically 
studied the effect of a first order chemical reaction on steady 
state dispersion of air pollutants from a point or infinite line 
source under atmospheric inversion. Nunge and VaidyanathanCl977) 
have also obtained analytical solution to the time dependent 
advective diffusion equation for dispersion of a reactive 
pollutant. Naresh and Nath(1990) have studied the effect of 
variability of chemical reaction on a two dimensional steady state 
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dispersion of air pollutant from ground based area sources with 
variable wind and diffusivity profiles and concluded that the 
concentration of the pollutant decreases as the rate of chemical 
reaction increases. 

1.5.4 Combined effects of chemical reaction, deposition and 
gravitational settling 

The studies related to deposition of pollutants on the ground 
has also been taken up. Since it is expressed as a flux, it can be 
readily incorporated into a solution of advective diffusion 
equation as one of the required boundary condition. But in the 
case of Gaussian plume model, since it does not explicitly treat 
the vertical transport process, it is difficult to modify it for 
deposition that is realistic. Removal of pollutant by deposition 
in the Gaussian model was discussed by Chamber lainC 1953 ) as noted 
in Van der Hoven(1968) referred to as the source-depletion model. 
Chamber lain( 1953) assumed that the shape of the plume is not 
affected by the deposition but only its total content. The 
limitations of this method has been discussed by many authors 
e.g., HorstC1977), Prahm and BerkowiczC 1978) , Nordlund et 
a/. (1979). The source-depletion model does not treat the vertical 
diffusive transport in accordance with the gradient transfer 
theory and neglects the effects of variation of diffusivity and 
wind with height. Assuming linear plume diffusion with downwind 
distance, Jensen(1980) presented relatively simple model including 
combined effect of diffusion and deposition. His, so called 
source-depletion model, however, becomes inferior to Horst ’ s C 1977 ) 
surface-depletion model especially under conditions where the 
atmospheric diffusion is slow and the deposition is fast. With 
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these limitations he gave explicit formulae for downwind 
deposition amounts and ground concentrations. CorbettC 1981 ) 
discussed the source-depletion method, which is used to account 
for dry deposition in Gaussian plume models of the atmospheric 
dispersion of pollutants, on the basis of a diffusivity model and 
used these results to derive simple criteria for the validity of 
the method, which can be applied within the framework of the 
Gaussian model without explicit recourse to diffusivity theory. 
Bartz i s ( 1 983) extended the Gaussian model by incorporating the 
deposition effects for dispersion of radioactive material released 
from an elevated point source. His model, based on ideal diffusion 
problem, is rather more realistic than the widely used 
source-depletion model. 

Numerous studies have been conducted to properly understand 
the phenomenon of deposition both mathematically and 
experimentally by taking into account the effects of deposition 
and also chemical reaction. Scriven and FisherC 1975a, b) , Rasmussen 
et ai.(1975), HogstromC 1975) , Mcmahon et ai. (1976), 
OvercampC1976) , Slinn(1976), Ermak(1977), Gar 1 and( 1977 , 1978 ) , 
Horst C 1977 , 1984) , Berkowicz and Prahm(1978), Platt(1978), Astarita 
et a/. (1979), Sehmel ( 1980) , Alam and Seinf eld(1981 ) , Davidson et 
ai .( 1981 , 1985) , Carmichael and Peters(1984a,b) , Renner et 
ai. (1985), Bennet( 1988) , Shukla and- Chauhan( 1^88 ) , Ch:auhB.n and 
Shukla(1989) , Koch(1989), Lehmann( 1989) , Noll and Fang(1989), 
Lovett and Kinsman(1990) , Pandis and Seinf eld( 1990) , 
Underwood( 1991a, b) , Rao et ai.(1992). The effect of gravitational 
settling of pollutants has also been studied e.g., Calder( 1961 ) , 
Smith(1962), Cleary et ai.(1974), Ermak(1977), Peterson and 
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Seinfeld( 1977) , Fisher and MacqueenC 1 98 1 ) , Rao ( 1 982 , 1 984 ) , 
Lee(1985), Rao and Godowi tch( 1 990 ) . The dispersion of air 
pollutants when ground absorption and gravitational settling can 
not be ignored has been discussed by Calder ( 1961 ) . In his 
formulation both the gravitational settling flux and the ground 
deposition flux are taken to be proportional to the local air 
concentration. Using this type of formulation, solutions to the 
atmospheric transport equation have been obtained for a number of 
cases. RoundsC 1955) , Smith(1962), Heines and PetersC 1974) . Using 
constant diffusivity and wind, Monin(1959) obtained an analytical 
solution of the diffusion equation including dry deposition. 
Scriven and FisherC1975) gave an analytical solution for the 
cross-wind integrated concentration from an elevated continuous 
point source taking into account the dry deposition of the gaseous 
pollutant at the surface and its washout in rain. ErmakCl977) has 
solved the steady state three dimensional transport equation for 
an elevated continuous point source and obtained analytical 
solutions incorporating dry deposition and gravitational settling 
of the pollutants. Fisher and MacqueenC 1981 ) have obtained steady 
state distribution for the integrated cross wind concentration of 
dust particulates emitted from an elevated point source including 
deposition and gravitational settling under inversion layer. 
Koch(1989) gave an analytical - solution of the steady state two 
dimensional atmospheric diffusion equation with dry deposition 
using power law profiles for wind velocity and vertical diffusion 
coefficient in an unbounded atmosphere while LehmannC 1989) 
presented a solution of time dependent one-dimensional diffusion 
equation for an inert pollutant from an elevated point source 
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incorporating dry deposition at the ground and leakage at the top 
of the inversion layer for simple diffusion coefficient profiles 
and illustrated them numerically. Rao and GodowitchC 1990) have 
briefly described a methodology to include the combined effect of 
deposition, gravitational settling and first order chemical 
transformation in urban air pollution models. The concentration 
expressions given by him can be thought of as analytical 
extensions of the familiar Gaussian plume dispersion models with 
these removal mechanisms. 

It is well known that a primary pollutant is often converted 
to a secondary pollutant in the atmosphere, (e.g.,S 02 is converted 
to SO^ ) which is even more deleterious. In this direction, Alam 
and Seinf eldC 1981 ) have studied the dispersion of sulfur dioxide 
and sulfate from a point source by taking into account the dry and 
wet removal of both the species and obtained analytical solutions 
of the steady state, coupled, three dimensional atmospheric 
diffusion equation incorporating first order chemical conversion 
and removal under inversion conditions. Similar problem including 
combined effect of settling, deposition, first order removal and 
conversion of gaseous to particulate pollutant has been studied by 
Peterson and Seinfeld (1977). Lee(1985) also gave analytical 
solutions for both primary and secondary pollutants including 
above mentioned factors in an unbounded atmosphere i.e., without 
inversion. 

It is pointed out here that these problems , where deposition, 
sedimentation, etc. are taken into account, have been studied 
under steady state conditions only but the sources of the 
pollutant could be time dependent e.g., explosion of a nuclear 
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device, leakage of toxic gases from storage tanks, etc. In such 
cases unsteady state atmospheric diffusion equation with above 
removal mechanisms should be considered, e.g., Cleary et 
a7.(1974), LI ewelynC 1983) , Shukla and ChauhanC 1988) . In 
particular, Shukla and ChauhanC 1988) have solved unsteady state, 
coupled, three dimensional convective diffusion equation for a 
reactive pollutant emitted from a time dependent point source 
forming secondary species in the atmosphere by taking into account 
the dry and wet deposition of pollutant for different forms of 
flux at the point source and applied the analysis to study the 
dispersion of MIC leakage at Bhopal, India, 

A thorough review of long range transport models and 
processes with regard to dry deposition and other removal 
mechanisms have also been described by some workers which are of 
special interest may be noted, e.g., Scriven and FisherC1975a,b) , 
FisherC 1975 , 1981 , 1983) , Eliassen (1980), Fay and RosenzweigC 1980) , 
SI innC 1982b) , Fisher and Smith(1987), Janssen and Asman(1988). 

In view of the above literature survey, it may be pointed out 
that the unsteady state dispersion of pollutant in the atmosphere 
has not been studied in a patchy atmosphere. Therefore, in 
chapters II and III of the thesis, we have studied the unsteady 
state dispersion of a reactive air pollutant from a time dependent 
point source forming secondary species for different forms of flux 
at the source by taking into account the first order chemical 
conversion of primary pollutant to secondary pollutant, removal 
process, deposition at the ground, leakage at the inversion 
height, gravitational settling, etc. in a two patch environment 
under inversion conditions. 
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1.5.5 Effect of greenbelt 

The concentration of pollutants in the atmosphere can be 
reduced considerably by plantation of greenbelt suitably at a 
required place to protect the region under consideration (a 
historical monument, a habitat). It is well established that 
certain plants have the capability of removing the pollutants by 
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that the species planted be relatively resistant to the exposure 
of particular pollutant. Several studies to study the effect of 
various pollutants on different plant species have been made in 
the last few decades or so, e.g., Thomas and Hendr icks ( 1956) , Todd 
and Garber ( 1958) , Brandt and HeckC1968), TreshowC 1968) , 
DainesC 1968) , Davis(1972), Garsed et a/. (1981), Takemoto and 
Noble(1982), Dochinger ( 1 988 ) , Evans(1988). It may be noted that 
there are plants which are relatively resistant to specific 
pollutants. In particular, survey of Indian plants in relation to 
atmospheric pollutants have been conducted by Ahmad et a/. (1988). 
They have found forty eight species of pollutant tolerant trees, 
shrubs and herbs that thrive well in the polluted environments and 
have been recommended for landscaping urban and industrial areas 
for mitigating air pollution. Some of them are listed here, 
Neem(Azad 2 rac/jta indica') , TLadambi^Anthocephalus cadamba) , Shisham 
(Dalbergia sissoo) , TenduCDiospyros melanoxylon') , Sahjan(Moring‘a 
oleifera') , AshokCPolyal thia 1 ongifol ia) , JamnniSyzygium cuminii), 
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K&ner (.Nerium odorum) , Amaranthus graecizans. 

The transport of gases and particulate matters within plant 
and vegetative canopies has been studied by many researchers, 
MacBean(1968) , Hill(1971). Bennet and Hi 1 1 C 1973 , 1 975 ) . Hartog and 
Shaw(1974), Raynor et a7.(1974), Shaw et ai.(1974), Aylor(1975). 
Chamberlain(1975a,b) , Petit et aJ.(1976), O’Dell et a7.(1977), 
Bache( 1979a, b) , Smith(1981), Hosker and Lindberg( 1982) , 
SI inn( 1982a) , Kapoor and Gupta(1984), Wiman and Agren(1985), Wiman 
and Lannef or s ( 1 985 ) , Wiman et a/. (1985), Baldocchi et aJ.(1987), 
Shukla and ChauhanC 1987 ) , Lindberg and LovettC 1992) , Lovett and 
LindbergC 1992) . In particular, Petit et ai.(1976) presented some 
results concerning characteristics of air flow within and above a 
forest and calculated SO 2 fluxes at the top of the canopies. He 
aimed at determining the behavior of the forest towards SO 2 and to 
develop a diffusion model of SO 2 with associated absorption by the 
roughness. Bache(1979a,b) has also developed an analysis and 
suggested a modified form of turbulent diffusion equation to take 
into account the particulate transport behavior within plant 
canopies. A theoretical framework to predict particle deposition 
due to vegetation by considering variable wind velocity profile 
has been given by SI inn( 1982a) . A review of atmospheric deposition 
regarding plant assimilation of gases and particles has been 
presented by Hosker and LindbergC 1982) and Smith(1981 ). Wiman and 
Agren(1985) proposed a mathematical model for aerosol depletion 
and deposition on forests by considering a modified form of 
convective diffusion equation. His model considers the interaction 
between forest structure, open field, forest aerodynamics and 


aerosol characteristics. 
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From these studies it may be speculated that if a greenbelt 
(in the form of forest capable of absorbing the pollutant) is 
provided suitably between the source and the region to be 
protected, it is possible to minimize the effect of pollution on 
the region under consideration. In this direction Kapoor and 
GuptaCl984) have studied the attenuation of an inert pollutant by 
a greenbelt under steady state conditions. 

It may be noted here that no information is available on 
concentration reduction due to greenbelt plantation when the 
pollutant is emitted from time dependent sources as well as the 
effect of inter-exchange process near the canopy surface. In view 
of this in chapter IV, the effect of greenbelt on the 
concentration reduction of a reactive air pollutant emitted from a 
time dependent point source with different forms of flux is 
studied. Also the optimum location of greenbelt between the source 
and the receptor (i.e., the area to be protected) is found for the 
maximum concentration reduction at the receptor. 
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1.5,6 Effect of precipitation scavenging 

In the last few decades increased number of 
have also been carried out to understand the 
pollutants/contaminants by precipitation. Engelmann 
(1970), Engelmann and Slinn(1970), Hales 
Hales et a/. (1973), Dana et al . (1973, 1975) . 

Davies( 1976) , Peters( 1976) , Barr ie( 1978) , 
Maul (1978), Radke et a7.(1980), Rodhe(1980), 

et 
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Hales(1972) presented some fundamentals for the general analysis 
of precipitation scavenging emphasizing the importance of 
reversible phenomenon. Since the process of trace gas scavenging 
is reversible in nature and the phenomenon of absorption and 
desorption may cause a redistribution of pollutants in the 
atmosphere as pointed by Hales(1972), Slinn(1974) has analytically 
considered the redistribution of gas plume caused by reversible 
washout and presented solutions in some simple cases. A simplified 
linearised model for predicting reversible gas washout of 
pollutant from industrial plumes has been presented by Hales et 
a7.(1973) considering rain to be composed of spherical, 
non-interacting drops of static size distribution. Peters(1976) 
presented a theoretical analysis of the contribution to rain water 
sulfate concentration by precipitation scavenging of gaseous SO 2 
and sulfate containing aerosols. The reversible rain washout of 
soluble gases such as sulfur dioxide, Ammonia and carbon dioxide 
has been studied by a number of authors, e.g., Dana et 37.(1975), 
Hill and Adamowicz C 1 977 ) , Adeuwyi and Carmichael ( 1982) . Shukla et 
ai.C1982) have calculated vertical concentration distribution of a 
gaseous pollutant emitted at ground level into a stagnant fog 
layer containing uniform distribution of fog droplets by assuming 
reversible absorption of a pollutant species. Fisher(1982) has 
studied the transport and removal of sulfur dioxide in rain using 
a simple model which combines the microphysics of absorption and 
chemical transformation of sulfur dioxide in clouds with the 
dynamics of the air motion and obtained approximate one 
dimensional solution of the governing equations. A detailed theory 
of HNOg vapor removal (particularly, below-cloud removal by rain) 
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has been carried out by Levine and Schwar tz ( 1 982 ) . Chang(1984) has 
derived wet removal coefficients for HNO^ vapor in rain and snow 
systems and parameterized them in terms of precipitation rate, 
under a number of approximations. Since the concentration of the 
pollutant in gaseous phase directly influence the concentration of 
the pollutant in droplet phase it is, therefore, necessary to 
consider simultaneously the coupled process of gas phase depletion 
and aqueous phase accumulation of the pollutant. Kumar(1985) has 
given an Eulerian model to describe the simultaneous process of 
trace gas removal from the atmosphere and absorption of these 
gases in rain drops by considering the precipitation scavenging of 
these gases present below the cloud base but his model does not 
take into account the wind velocity and the presence of a 
continuous source. He (1986) has extended the above model to take 
into account the process of absorption of multiple species and 
chemical reactions within the droplets. 

In the above studies the transient emission of pollutant from 
the source as well as the non-linear interaction between the 
primary and secondary pollutants with the time dependent 
distribution of rain droplets have not been taken into account. In 
view of this in chapters V-VI , we have studied the effect of 
precipitation scavenging on concentration distribution of a 
reactive air pollutant from a time dependent point source for 
different forms of flux at the source. In these studies we have 
considered two phases in the atmosphere i.e., the pollutant phase 
and the phase of absorbed pollutant in the droplets whereas in 
chapter VII we have proposed a nonlinear ecological type model to 
study the effect of precipitation scavenging on the concentration 
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distribution of primary and secondary pollutant using stability 
theory . 

1.6 SUMMARY OF THE WORK DONE IN THE THESIS 

The thesis consists of seven chapters and deals with the 
unsteady state dispersion of a reactive air pollutant emitted from 
a time dependent point source with different forms of flux at the 
source. The main concern has been to discuss the effects of 
removal mechanisms such as chemical transformation, dry 
deposition, removal by greenbelt and precipitation scavenging on 
the dispersion process. In these studies, the unsteady state 
convective diffusion equations both for primary and secondary 
pollutants have been solved with appropriate initial and boundary 
conditions. The effect of patchiness in the atmosphere on the 
dispersion process has also been studied. It is assumed that the 
wind is sufficiently large to neglect the downwind diffusion. A 
nonlinear model is also proposed for removal mechanism of both 
primary and secondary pollutants due to precipitation scavenging 
using ecological type ideas in the last chapter of the thesis. 

In chapter I a general introduction and basic concepts of air 
pollution with relevant literature is presented to provide a 
necessary feedback required for the forthcoming chapters. 

In chapter II the unsteady state dispersion of a reactive air 
pollutant emitted from a time dependent point source forming 
secondary pollutant is studied by taking into account the first 
order conversion of primary pollutant to secondary pollutant in a 
two patch environment, assuming the existence of sinks, uniformly 
distributed throughout the second patch, causing the removal of 
both primary and secondary pollutants. These sinks in the 
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here may represent rain droplets, fog droplets or even 
ally introduced species (liquid or gas) which are effective 
removal of pollutants by chemical reaction, etc. The latter 
s may be useful in removing the toxic species in the 
lere by accidental releases from an instantaneous source, 
cblem is studied by solving unsteady state, coupled, three 
onal convective diffusion equations in both the patches for 
Drimary and secondary pollutants with removal terms 

rating the deposition at the ground and the leakage at the 
on layer with suitable matching at the interface of the two 
. The time dependent flux at the source is considered to be 
forms, Ci) Instantaneous (u) Constant, and (iii) Step 
1 type. 

5 effects of removal mechanisms on the concentration 
ition of both the species are then discussed in all the 
ises of time dependent flux at the source. It is observed 
le concentration of pollutants decreases as downwind 
increases in all the three cases. Also the concentration 
iry pollutant decreases but that of secondary pollutant 
s at a particular time and location as the rate of 
on of primary pollutant to secondary pollutant increases, 
ound that the central line concentration of both the 
ts decreases as time increases for instantaneous flux 
it increases as time increases, and reaches their 
re steady state values for constant flux. When flux is 
;tion type, the concentration characteristics of both the 
s are similar to the combined effect of the above two 
is noted that the concentration of both the pollutants 
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decreases considerably due to removal parameter associated with 
the sink phase. 

From the analysis it is concluded that the presence of sink 
phase in the second patch in the form of rain droplets, fog 
droplets or even externally introduced species can be very 
effective in removing the pollutants from the atmosphere. The 
analysis thus suggests a mechanism by which toxic gases leaked out 
in the atmosphere due to accidental discharge, etc. can be 
neutralized by introducing suitable liquid or gas phase in the 
atmosphere by mechanical means. The analysis can also be utilized 
for studies related to risk analysis particularly in cases of 
accidental leakage of toxicants/pollutants such as MIC leakage in 
a pesticide plant, Bhopal, India. 

In chapter III we study the unsteady state dispersion of both 
primary and secondary pollutants, which are heavier than air, in a 
two patch environment to see the effect of gravitational settling, 
on the dispersion of these pollutants. The deposition of settling 
particles at the ground and the removal of these pollutants have 
been taken into account. As in chapter II, the relevant equations 
with settling terms have been solved and it is found that the 
concentration of the pollutants decreases when gravitational 
settling velocity increases for the points above the source 
height, however, for the points below the source height, result is 
just reversed. This is expected because after emission of 
pollutants from the source gravity comes into effect and heavier 
particles are likely to have higher settling velocity in 
comparison to not so heavier particles. The effect of removal is, 
however, to decrease the concentration of both the pollutants. 
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Chapter IV deals with the effect of greenbelt (plantation 
capable of absorbing the pollutant and inducing deposition) on 
dispersion of a reactive air pollutant from a time dependent point 
source under atmospheric inversion condition. The inversion layer 
is divided into four regions where source is located in the first 
region and the greenbelt lies some distance away from the source 
in downwind direction in the second region. The region above the 
greenbelt is the third region whereas the region next to greenbelt 
in downwind direction is the fourth region i.e., the region, where 
the receptor (a historical monument, a habitat, etc.) is located. 
The vertical and horizontal diffusion processes in the greenbelt 
region are assumed to be negligible due to denseness of the plant 
canopies and shrubs grown on the ground below canopies of the 
plants. It is assumed further that the entire greenbelt is acting 
as a sink for the pollutant under consideration. The wind 
velocity, diffusion coefficients and deposition rates are taken to 
be different in all the four regions. The relevant equations with 
above assumptions have been solved within and outside the 
greenbelt using appropriate boundary and matching conditions with 
suitable deposition condition at the top of the greenbelt. The 
analysis has also been done for the case where inter-exchange 
process is taken into account, in a narrow zone of the third 
region adjacent to the top of the greenbelt, characterized by an 
exchange velocity. It is found that the concentration of the 
pollutant decreases considerably in presence of greenbelt. The 
optimization of the location of greenbelt is also studied to 
affect the maximum concentration reduction on the receptor in the 
above cases. It is seen that, when the source height is greater 
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han the greenbelt height, the greenbelt should be located close 
o the receptor and when the source height is less than or equal 
o the greenbelt height, it should be located near the source 
•e., if sources are ground based like storage tanks it is 
esirable to surround this area by the greenbelt to avoid the 
ollutant spreading in the case of an accidental blast. The effect 
f other removal parameters like deposition, depletion coefficient 
ue to greenbelt, inter-exchange velocity, etc. is to decrease the 
oncentration of the pollutant. 

Chapter V deals with the study of the effect of precipitation 
cavenging on the dispersion of a reactive pollutant emitted from 
time dependent point source. It is assumed that the rain 
roplets are uniformly distributed in the atmosphere which absorb 
he pollutant and remove it by their fall on the ground with a 
onstant velocity. We consider that the concentration of the 
ollutant in the gaseous phase is governed by the unsteady state 
hree dimensional convective diffusion equation with removal terms 
hereas the concentration of the absorbed pollutant in droplet 
base is assumed to be governed by a simplified equation which 
akes into account the falling of rain droplets but neglects the 
if fusion, the case applicable to the heavy monsoon rain. These 
quations are then solved using suitable initial and boundary 
onditions and the expressions are obtained for the concentration 
istribution of pollutant in gaseous phase and in rain droplet 
base. The forms of time dependent flux at the point source are 
aken to be Ci) Instantaneous, and (ii) Constant. It is shown that 
he concentration of the pollutant in the atmosphere reduces 
onsiderably, under appropriate conditions, due to precipitation 
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scavenging by rain, thus providing an important mechanism for 
removing pollutants/toxicants from the atmosphere. 

In chapter VI a similar problem of removal of pollutant by 
precipitation scavenging is studied in a two patch environment. It 
is assumed that the rain is taking place, some distance away from 
the source in the second patch where the effect of precipitation 
scavenging is studied. It is shown that, as the pollutant front 
reaches in the second patch where rain is taking place, the 
concentration of pollutant reduces significantly due to washout 
process for both the cases of instantaneous and constant flux at 
the source. 

Finally, from the above analysis and the results in chapter 
V, it may be concluded that precipitation scavenging is a very 
important mechanism for the removal of pollutant from the 
atmosphere at a given receptor point located in the second patch. 
However, the rate of reduction depends upon the location of the 
patch and prevailing meteorological conditions. It is also noted 
that the reduction in pollutant concentration is larger at this 
receptor when precipitation patch is closer to the source. 

Since in real situations during precipitation, the number 
density of rain, fog or snow droplets changes as the intensity of 
precipitation increases and this change in number density affects 
the interaction process and should, therefore, be taken into 
account in the model. Thus, In chapter VII a nonlinear model is 
proposed for removal mechanism of both primary and secondary 
pollutants, due to precipitation scavenging, by taking into 
account the interaction of the droplet phase with the pollutant 
phase. In fact, we consider that there exist four separate 
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nonlinearly interacting phases, (i.e., the droplet phase, primary 
pollutant phase, secondary pollutant phase and the combined phase 
of both the pollutants absorbed in droplets). We assume that in 
these phases, the number density of droplets, the concentration of 
both primary and secondary pollutants and the concentration of 
both the pollutants absorbed in the droplets are governed by 
evolution equations (i.e., nonlinear ordinary differential 
equations) with source and removal terms. It is assumed that the 
rate of removal of pollutants is proportional to the concentration 
of pollutants as well as the number density of droplets in the 
atmosphere. Further, the concentration of pollutants absorbed 
depends nonlinearly on the number density of droplets as well as 
the concentration of both the pollutants. Our aim here is to 
analyze the proposed nonlinear model and to see the effect of 
removal mechanism on the equilibrium level of pollutants in the 
atmosphere. The analysis is done for both cases of with and 
without secondary pollutant. To understand the removal mechanism, 
Liapunov’s stability is used. It is shown that when the pollutant 
is emitted by an instantaneous source (single introduction of the 
pollutant in the atmosphere), both the pollutants would be removed 
completely from the atmosphere by the precipitation process. When 
the pollutant is emitted at a constant rate, it is found that the 
pollutants may be removed from the atmosphere completely but it 
would depend upon the rate of emission and other removal 


processes . 



CHAPTER 11 


Unsteady State Dispersion of a Reactive Air Pollutant from a Time 
Dependent Point Source Forming Secondary Species in a two Patch 

Environment 

2 . 0 introduction 

In the environment various mechanisms, such as chemical 
reaction, dry deposition on the ground and on vegetative canopies 
and precipitation scavenging by rain, fog or snow, are known to 
exist and are effective in the removal of pollutants from the 
atmosphere. In the last few decades several investigations have 
been conducted to study the dispersion of air pollutants in the 
atmosphere under various meteorological conditions with different 
kinds of removal mechanisms, Scriven and Fisher(1975a,b) , McMahon 
et ai.(1976), Alam and Seinf eld( 1981 ) , Shukla et ai.(1982), 
Robson(1983) , Llewelyn(1983) , Seinfeld(1986) , Shukla and 
ChauhanC 1988) , Chauhan and Shukla( 1989) , Lehmann( 1989 ) , Shukla et 
ai. (1991b) including scavenging by rain droplets, HalesCl972), 
Hales et a). (1973), Slinn(1974), Kumar ( 1985) . 

It is known that primary pollutant is often converted to 
secondary pollutant in the atmosphere (e.g., SO 2 is converted to 
SO^ ) and both of these are required to be removed from the 
atmosphere. In this direction Alam and Seinf eld( 1981 ) have studied 

Part of this chapter is to appear in the Proceedings of Fifth 
International Conference on Precipitation Scavenging and 
Atmosphere-Surface Exchange Processes , 15-19 July ,i99l,Ri chi and , 
Washington, U.S.A. 


f y 


37 


the dispersion of sulfur dioxide and sulfate from a point source 
by taking into account the dry and wet deposition of both the 
species and obtained analytical solution of steady state, three 
dimensional atmospheric diffusion equation incorporating the first 
order conversion and removal, Shukla and ChauhanC 1 988 ) have also 
studied the similar problem under unsteady state conditions when 
the pollutant is emitted from a point source, with different forms 
of time dependent flux and applied the analysis to the study of 
dispersion of MIC leakage in a pesticide plant in Bhopal, 1984, 
India . 

As pointed out in chapter I, in real situations, the 
ecological conditions on the ground and meteorological conditions 
in the atmosphere are different at different places leading to 
patchiness in the environment with different dispersal , deposition 
and removal parameters in each patch. This aspect has not been 
studied in the above investigations. 

In this chapter we, therefore, study the unsteady state 
dispersion of air pollutant emitted from a time dependent point 
source forming secondary species by taking into account the 
effects of removal mechanism in a two patch environment, assuming 
the existence of sinks uniformly distributed throughout the second 
patch causing the removal of both the primary and secondary 
pollutants. In the atmosphere these sinks may represent rain 
droplets, fog droplets or even externally introduced species 
(liquid or gas) which are effective in the removal of pollutants 
by chemical reactions, etc. The latter process may be useful in 
removing the toxic species in the atmosphere by accidental 
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releases from instantaneous sources. This problem has been studied 
by solving coupled three dimensional unsteady state convective 
diffusion equations involving both primary and secondary 
pollutants with removal terms and incorporating the deposition at 
the ground and the leakage at the top of the inversion layer with 
suitable matching at the interface of the two patches. The forms 
of time dependent flux at the point source considered in the 
analysis are: (i) Instantaneous, (ii) Constant, and (iii) Step 
function type. 

The effects of removal mechanism in the second patch on the 
concentration distribution of both primary and secondary species 
are discussed in all the above cases. 

2 . 1 MATHEMATICAL MODEL 

Consider the dispersion of a reactive air pollutant, 
undergoing first order chemical conversion and removal, from a 
time dependent point source of strength QCt) located at height h 

s 

from the ground z =0, forming secondary species in a two patch 
environment under atmospheric inversion as shown in Figure 2.1. 
The source of secondary pollutant is, however, the transformation 
of primary pollutant. It may be noted that the time dependent 
source is located in the first patch where some kind of removal 
mechanism may exist, however, in the second patch it is assumed 
that there exist a phase of uniformly distributed sinks as 
described earlier which are mainly responsible for removing both 
the primary and secondary pollutants from the atmosphere. In the 
two patches wind velocity, deposition and other removal parameters 
are taken to be different. It is also assumed that the wind is 
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sufficiently large in both the patches so that downwind diffusion 
is negligible in comparison to advection. 

Under these assumptions the partial differential equations 
governing the concentration of pollutants in both the patches, 
with appropriate initial and boundary conditions can be written as 
f ol lows , 


Patch I (0£x<Xj, 0:SziH) 

The differential equation governing the concentration 
Cj^(x,y,z,t) of the primary pollutant in the first patch is written 


ac^ aCj a^Cj 

dt~ * "" ^yl TIT" 

ay 




( 2 , 1 ) 


where x, y, z are cartesian co-ordinates, t is the time, x-axis is 
taken in the downwind direction and z-axis vertically upwards, u^ 
is the mean wind velocity, and are dif fusivities in y and 

z directions respectively, taken as constants. The constant k^ is 
the conversion rate coefficient of primary pollutant to secondary 
pollutant and k^, is its removal rate coefficient. 

We prescribe the following initial and boundary conditions 
for equation (2.1), 

The initial concentration distribution of the pollutant is 


assumed to be zero. 


Cj(x,y,z,t) =0 at t = 0 (2,2) 

A point source of strength Q(t) is assumed to be located at 

(0,0, h ). Thus the condition at x = 0 is, 
s 
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C.(x,y,z,t) = ^ ^ 5(y) 6(z-h ) at x = 0 (2.3) 

1 t» s 

At infinite lateral distance the concentration approaches 

zero , 

Cj(x,y,z,t) =0 as y — > ± oo (2.4) 

The pollutant may also be removed by the deposition process 
at the ground z = 0, 

ac 

hi W- ‘ ‘It ^ = 0 (2.5) 

The pollutant may also ‘escape’ into the base of the 
inversion layer at z = H, which is taken into account as, 

ac. 

^ at z = H (2.6) 

zi az 11 


where v^^ is the deposition velocity of the primary pollutant at 
z = 0, which in general depends on the particular species or 
material that is being removed and on the prevailing 
meteorological conditions and w^^ is its leakage velocity at z = H, 
the height of the inversion layer, a layer of stable air which 
forms a barrier to vertical mixing and 6( . ) is the Dirac delta 
function. 

Similarly the differential equation governing the 
concentration (x ,y , z , t ) of the secondary pollutant can be 
written as. 


dC 


pi 


ac 


at 


+ u 


1 ax 


PA = K 


yi 


+ F 

dy 


a^c 

— 7 ^ + k.c. - k^.c , 

3z 2 11 pi pi 


(2.7) 


The initial and boundary conditions for equation (2.7) are, 
(assuming that there is no direct emission of secondary pollutant 
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from the source) 


Cpi(x,y,z,t) =0 at t = 0 


C .(x,y,z,t) =0 at x=0;y=0, z=h 
X s 


Cpi(x,y,z,t) = 0 as y — > ± on 


dC 


K 


pi _ 


z 1 dz 
dC 


K 


pi _ 


zl dz 


V , , C , at z = 0 
dpi pi 


w . C . at z = H 
pi pi 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


where k , is the removal rate coefficient 
pi 

pollutant, v^p^ is its deposition velocity at z 
leakage velocity at z = H. 


of 
= 0 


the secondary 
and w , is the 


Patch II (x i x^, 0 z H) 

The differential equation governing the concentration 
C 2 (x,y,z,t) of the primary pollutant in the second patch is 
written as, 


ac ac a^c a^c 

* “2 35 r = ’‘ y 2 - «2 " ‘‘*2 " « 1>'=2 

oy dz 


(2.13) 


where U 2 is the mean wind velocity, Ky 2 dif fusivities 
in y and z directions respectively, k 2 is the conversion rate 
coefficient of primary pollutant to secondary pollutant, 1® 
its removal rate coefficient and 6^ is the depletion rate 
coefficient due to uniformly distributed sinks, the magnitude of 
which may increase if number or size of sinks increases. 

The initial and boundary conditions for equation (2.13) are. 


C2(x,y,z,t) = 


0 at t = 0 


(2.14) 
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C 2 (x,y,z,t) = Cj(x,y,z,t) at x = C2.15) 

C 2 (x,y,z,t )=0 asy — >±oo (2.16) 


ac2 

^z2 ^ ''d2^2 


at z = 0 


(2.17) 


K 


ac2 

z2 az 


at z = H 


(2.18) 


where v ^2 is the deposition velocity of the primary pollutant at 
z = 0 in the second patch and W 2 is its leakage velocity at z = H. 

Similarly the differential equation governing the 
concentration Cp 2 (x,y,z,t) of the secondary pollutant can be 
written as, 


ac^„ ac^„ 

P2 . ,, p2 _ ^ 

at 2 ax ■ 


y2 ^ " ^Z2 ^ " “2^ - 'V2"«2>Cp2 '2.19) 


The initial and boundary conditions for equation (2.19) are, 


Cp 2 (x,y,z,t) =0 at 

0 

II 




(2.20) 

CpjCx.y.z, 

t ) = . ( X 

pi 

.y.z, t) 

at 

X = Xj 


(2.21) 

Cp2(ii.y.2. 

t ) = 0 as 

y — » ± 

00 



(2 . 22) 

sc„, 

XT p2 ^ 

z2 az 

''dp2S2 

at z = 

0 



(2.23) 

XT p2 _ 

z2 az 

'^p2^p2 

at z = 

= H 



(2 . 24) 

where k ^ 

p2 

is the 

removal 

rate 

coefficient of 

the 

secondary 

pol lutant , 

§2 is the depletion 

rate coefficient 

of 

secondary 

pol lutant 

due to sinks, v, » 

dp2 

is the deposition velocity 

of it at 
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z = 0 and w „ is its leakage velocity at z = H. The following 
forms of flux Q(t) at the source are considered in the subsequent 
analysis , 

(i) Flux is Instantaneous 

Q(t) = 5(t) 

(ii) Flux is Constant (2.25) 

Q(t) = Q (constant) 
c 

(iii) Flux is Step function type 

Q(t)=Q^ O^tit 
c o 

= 0 t > t 

0 

where Q is constant. 

0 

It is convenient to cast the problem in dimensionless form. 
For that we use the following dimensionless quantities, 


K t 

T = 


H 


2 ’ 


X = 


K X 
zo 

u 

0 


— y — z T- 
S' = ^ " -ff-' ’'s = 


h u . 

® TT - ^ 

'TT ’ - — 


u u „ K 

r r t R c o(t)~^^ ^ o — Q ^ ® 

i -57 '^i- '"pi Q, Sr Q,, ' S ^2- 


hS. 

zo 


i = 1 ,2 


where K and u are reference diffusion coefficient and wind 
zo 0 

velocity respectively. 

The equations (2.1)-(2.25) can be written in dimensionless 
form (dropping bars for convenience) as : 


Patch I (0 ^ X < Xj , 0 i z i 1) 




= 


dy 


2 - * H 


a c. 


az' 


'“1 * “gi’°i 


aCj aCj 

at * ax 


(2.26) 
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Cj(x,y,z,t) = 0 

at t = 0 


(2.27) 

C,(x,y,z,t) = ^ 5(y) 5(z-h ) at 

Jl ^ J. ^ 

X = 0 

(2.28) 

Cj(x,y,z,t) = 0 

as y — > ± <0 


(2.29) 

ac. 

az = ''iS 

at z = 0 


(2.30) 

ac, 

aJ^ = - 

at z = 1 


(2.31) 

„ "Si _ 

at 1 ax 

a^c a^c 

S V " >1 ? * 

ay ^ az^ 

«1^1 " «pl^pl 

(2.32) 

Cpj (x,y,z,t) = 0 

at t = 0 


(2.33) 

Cpj(x,y,z,t) = 0 

at x=0;y=0, z 

= ^ 

(2.34) 

Cpj (x ,y , z , t) = 0 

as y — * ± CD 


(2.35) 

= N C 

az pi pi 

at z = 0 


(2.36) 

pA - _ t c 
az pi pi 

at z = 1 


(2.37) 

Patch II (x x^ 

, 0 :S Z ;£ 1) 



ac ac 

at ■*■ '^2 ax ^2 

a% 

7 ^ ’'2 ■ 2 

ay az 

+ «g2 + 5 i)C2 

(2.38) 

C2(x,y,z,t) = 0 

at t = 0 


(2.39) 

C 2 (x,y,z,t) = Cj(x.y,z,t) at x = Xj 

(2.40) 

C2(x,y,z,t) = 0 

as y — > ± 00 


(2.41) 



45 


ac„ 

£ = N c 

dz 2^2 


at z = 0 


(2.42) 


ac. 


-t. __ 

az 

4^ 

^S2 , 

5 

at “ 

2 ax 


z , t ) 

S2'*'>'’ 

z , t ) 


z , t ) 


p2 


at z = 1 


Po — ■*■ H'' — ® 


2.2 

ay 


'2 az2 


2^2 


"pi 


at X = X, 


(2.43) 


^«P2^^2^S2 ^2.44) 


(2.45) 

(2.46) 

(2.47) 


P^ = N C 

az ^p2 p2 


at z = 0 


(2.48) 


ac 


^ = - L C 

^2 p2 


dz 

where 


k.r 

«i = -K^‘ 
zo 


at z = 1 


Pi = 


a 


k^.H 
gi 


gi 


^Pi 


(2.49) 


K 

'’i et: ’ 

ZO 

Zl 

w . H 


pi 

r~7 

i = 1. 2 

Zl 


(i) 

Q(t) = 6(t) 

(ii) 

Q(t) = 1 

(iii) 

Q(t) = 1 , 0 :S t 


=0, t > t 


K . 

K . 

k .h2 

yi 

Zl 


K ’ 

K 

pi K ’ 

ZO 

zo 

zo 

.H 

V . .H 

w .H 

di 

K = ‘ipi 

T. = -i - 


pi 


z 1 


Zl 


(2.50) 
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2.2 METHOD OF SOLUTION 

Following Astarita et ai.(1979), Alain and Se inf e 1 d ( 1 98 1 ) , 
equations (2.26) and (2.32) can be written in compact form as, 




If we define a new concentration such that 



the uncoupled system for patch I can be written as, 


0 


(2.52) 
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- «pl®l = 0 


(2.53) 


oc. 


where B, = C . - 

1 pi Opj - aj - a 


C2.54) 


gl 


The conditions for remains same 

conditions for are. 


The corresponding 


Bj(x,y,z,t) =0 at t = 0 


(2.55) 


a, 


Bj(x,y,z,t) = 


Q(t) 


“pi ■ “i ■ “gl “i 


6(y) 6(z-h ) at X = 0 (2.56) 

s 


Bj(x,y,z,t) =0 as y — > ± oo 


(2.57) 


dB. - N. 

1 = N ,B, + Ei 1 

3z pi 1 a 


pi “l “gl 


“I'^i 


at z = 0 


(2.58) 


aB^ 

dZ 


- L . B. + — 
pi 1 oc 


^1 - "pi 


pi “1 " “gl 


“i^i 


at z = 1 


(2.59) 


Similarly we also get uncoupled system for Patch-II 
corresponding to equations (2.38) and (2.44), 


^2^2 ^“2 “g2 ^1 ^^2 ■ ® 


(2.60) 


^2®2 ^®p2 ^2^®2 ~ ° 


(2.61) 


dt 


where £0 ■ - ^ ^ + y 


i- + B ^ + X ^ 

2 ax ^2 ay2 ^2 az2 


(2.62) 


a. 


and B„ = C „ - ^ C„ 

2 p2 ap2 «2 ®g2 ^2 ^ ^ 


(2.63) 


The conditions for C 2 remains same, 
conditions for B 2 are. 


The corresponding 


B2(x,y,z,t) 


0 


at t = 0 


(2.64) 
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B2(x , y , z , t ) 


B2(x , y , z , t ) 


Bj(x,y,z,t) + 


Of, 




^gi 




at X 


= 0 as y — > ± 00 


Xj (2.65) 

( 2 . 66 ) 


5B, 

I 

dT 


= N ^B^ + 


%2 - ^2 


p2 2 ap2 - «2 “ “g2 ^ ^2 ■ ^ ^ 


a^Co at z = 0 


(2.67) 


aB. 

dT 


L „B„ + — 
p2 2 a_ 


L, - L „ 

^ ^ ^ at z = 1 

- + Oo “ o, 2 2 


( 2 . 68 ) 


p2 "2 “g2 "2 "1 

2.3 CONCENTRATION DISTRIBUTIONS 

The solutions for C^ , B^ and C 2 , B 2 are obtained using 
Laplace and Fourier transform techniques, (Carslaw and Jaeger, 
1959) in each cases as. 


2.3.1 CASE I WHEN FLUX IS INSTANTANEOUS AT THE SOURCE 


Q(t) = 6(t) 


Q, 


B^(x,y,z,t) = — 


C2(x,y,z,t) = — 


Q(x,y,z) 6 


M(x,y ,z) 6 


ft - 



r 

“ij 


ft - 

1 






“i 



(2.69) 


(2.70) 


(2.71) 


B2(x,y,2,t) 


Q, 


u. 


N(x ,y , z) 6 


ft - 

f!i 

^ 11 

1 

>■ “1 



(2.72) 


where 


P(x,y,z) 


2 ^1 

-y2/(4^,,) 

e 1 


e.p[- -ilp^F^Cz, 


F (h_) 


4n 




n=l 


u. 


X 
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Q(x, y,z) = 


2 ^1 
-y^/(4^X) 

e ^ «1 

I 5 — L 

, n=l * S' 

— X 


“pi**!? 


' " 1x1 " 

^ J'‘J -F 


H_(z) 


»1 2 

IT- -V^ 111 Il2 - "n'h.' 

1 m 


M(x ,y , z) = 


^-y"/4^(x) ^ ^ , 

> exp - ■ 

J4nB(x> ^ y '■ 


- I — ^-u:; J'" - ’‘iM -Q- 




F (h ) 1 

m s r 


r F (z)G (z)dz 
J m n 


N(x , y , z ) = 


-y /4pCx) » (■ ,a„o+5,+y^y 


mecx. 


— in - "dj -F 


J_(2) 




2 m s 
r P 


I 2 I I 22 - 


(ajj-aj2^ exp^ 




r s 


J F^Cz)J^(z)dz * ^ ^ ^2 ^ 

0 r=l ^ m=l 


(-r-i 


“gi * iiVi 


F (h ) 1 

-f_5- J F„(z)G^<z)dz 


F (z) = cos X Z + r — Sin X z 
n n X_ n 

n 
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^2 

G (z) = cos u z + — sin u z 
n (j 

'^n 


N 


H^(z) = cos P z + sin ^^z 

n ^n F ^n 

^n 


N 

J (z) = cos p z + ■ sin p z 
n n f n 

n 


2 

P = (cos X z + T — sin X z) dz 
n J n X n 

0 


N 


r 2 2 

= (cos u z + — sin u z) dz 
n J ^n ^n 


0 


pi N 

P . = 1 (cos F z + sin F z) 

nl J ^n ^ "^n 


"dz 


N 


r 1)2 2 

P_o = (cos pz + -£— sin p_z)^dz 
n2 J n n 

^ n 


Here X^, and are 

transcendental equations, 


the roots of the 


tan X = 
n 


tan = 


X„(L, + N, ) 
n 1 1 


tan 5 


n 


tan = 
n 



L,N, 

n 

1 1 

/In (1-2 

" N^) 

2 

^n " 

L 2 N 2 


* Npi> 

- 

^n 

"pi^i 

y (L o 
n p2 

" Np2> 

2 

Lp2 Np2 


f ol lowing 
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111 = J 


1,2 = j exp[- [ 


«1 «gl ~ %1 ~ 

_ — — J 


Hi - 1 <"2-" " 


I 22 = J 


2 _ .2 

r * ”gl ^1 I'J x'jdx' 


; (ag^x - a,,,)J„(z). 


I 24 = J 


«2 * «g2 ^1 * HH ~ HH ' “p2 ' ^2- | X- 

I ■ , ' ' '' ' — 1 1 y 


^11 ■ «Dl ■ " "'gl 


'12 «p2 ^2 " “2 " “g2 1 


a - ^ CL, - Si)N,i F^(l) 


Lpi + ^1 - ^pi>Npl 


" 2 n ^Pl - ^1 


^(x) 


^1 ^2 


ixi.^Cx-x^) 


C •-. ^ '■'<■•’ "L 

i i, i „ KANPU« 

-— ~" ^ 1 1 7 i 
4ec.^ No. A.' • 


CNp2 - ^2 ~ ^ 2 '’ V 

' L „2 ■" “ " ‘'p 2 '"p 2 
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^6n _ ^p2 ^2 

^ ■ ^p2 

and 5(t-x) is the Dirac delta function defined as, (Carslaw and 
Jaeger, 1941) 

5(t-x) =0, t < X 

= X £ t :£ X + e 

= 0 , t > X + e 

2.3.2 CASE II WHEN FLUX IS CONSTANT AT THE SOURCE Q(t) = 1 
Cj(x,y, 2 ,t) = -i- PCx.y.z) H(t 

Bj(x,y,z,t) = — ^ Q(x,y,z) H^t 

C 2 (x,y,z,t) = — ^ M(x,y,z) H^t 

B 2 (x,y,z,t) = — ^ N(x,y,z) H^t 

where H(t-x) is the Heaviside 
Jaeger, 1941) 

HCt-x) =0, t i X 

= — X < t < X + £ 
x+e 

= 1 , t i X + e 

From equations (2.73) - (2.76) it is observed that if t — > oo 
the expressions reduce to the case of steady state. It may also be 
noted that when t < x/u^ , the concentration of pollutants is zero 
i.e., the pollutant front has not reached that point. Also from 
equations (2.73) and (2.74) it is noted that if and L, 



function defined 


(2 . 73) 

(2.74) 

(2.75) 

(2.76) 

, (Carslaw and 


are 
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zero i.e., if there is no leakage into the inversion base and t — > 
the solutions for and Bj reduce to the same form of steady 
state solutions as obtained by Alam and Seinf eld( 1981 ) . 

2.3.3 CASE III WHEN FLUX IS STEP FUNCTION TYPE AT THE SOURCE 
Q(t) = 1 , 0 i t i t^ 

=0, t > t^ 


B. 


1 

P(x,y.z) 

H 

ft - -^1 

- H 

t — f — . 

"" 1 

H(t-to)j 

^1 

JLX 

L Ujj 



1 

^1 

Q(x,y ,z) 

Hj 


- Hj 

ft-t -- 
L 0 

^ 1 

H(t-t^)j 


C^Cx.y.z.t) = M(x,y,z)[H[t “ ( ^ ^ ]] 

- - ‘o - (-V7 * ^ 


(2.79) 


B^Cz.y.z.t) = 4- N(z,y,z)[H[t - (-II 4 

- - ‘o - ( ^ ^ 


(2.80) 


From equations (2 . 77 )-(2 . 80) it is seen that as t^ — > a>, 
these equations reduce to equations (2 , 73)-(2 . 76) respectively. 

The concentrations Cp^ and can now be obtained from C^ , 

B. and C^ , B„ in each of the above cases by the following 

1 jLl 


relations , 

Cp^(x,y,z,t) = Bj(x,y,z.t) + - 


a. 


pi " “l ~ “gl 


Cj(x,y,z,t) 
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Cp2(x ,y , z , t ) 


B_(x,y,z,t) + 

2 «p2 " «2 


“2 

«g2 ^2 ■ 


C2(x , y , z , t ) 


It may be noted here that the solutions above are written in 
terms of integrals for brevity only although these can be 
evaluated in exact forms which are very lengthy to write. 


2.4 DISCUSSION AND RESULTS 

To see the effects of removal parameters, the expressions for 

concentration distribution of both primary and secondary 

pollutants in both the patches have been computed and displayed 

graphically in Figures 2. 2-2. 8 . The parameter values used in the 

computation are taken as (Alam and Seinfeld, 1981): k^=0.05hr 

k^.=10 ^sec ^,k .=10 ^sec K .=5m^sec u . =2 . 25m sec \ 

gi ’pi ’ zi ’1 

w^= 0 . 1 cm sec Wp^= 0 . 01 cm sec \ (i = l, 2 ), v^ji=1.05cm sec 

V, 0 = 1' 0cm sec v, =0. 105cm sec v, o=0.10cm sec h =200m, 
dz dpi dp 2 s 

3 

H=10 m. The corresponding dimensionless values are found as 
a. =2. 77, ttg. =2.0. Op. =20.0, L.=0.2, Lp.=0.02, u.=1.0. Nj=2.1, 
N„=2.0, N =0.21, N „=0.20, h =0.20, 6 . =10.0, 7 . =1.0 and x =0.03, 
y=0.0, £=0.005. 

From the figures it is observed that the concentrations of 
both primary and secondary pollutants decrease as downwind 
distance increases in all the three cases. Also the concentration 
of primary pollutant decreases but that of secondary pollutant 
increases in all the cases at a particular time and location as 
the rates of conversion of primary pollutant to secondary 
pollutant i.e., and in the two patches increase. The effect 
of removal is, however, to decrease the concentrations of both the 
species in all the three cases. 
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When flux is instantaneous at the source, the central line 
concentrations (z = h , y = 0.0) of both primary and secondary 

S 

pollutants are plotted in Figures 2 . 2-2 . 3 for different values of 
t. It is noted that as time increases the concentration of both 
primary and secondary pollutants decreases. Comparing Figures 2.2 
and 2.3 it is seen that the concentration of both the pollutants 
reduced considerably as removal parameter, associated with the 
uniformly distributed sink phase, increases. 

When flux is constant at the source, the central line 
concentrations of both the pollutants are shown in Figures 2 . 4-2 . 5 
for different values of t. From these figures it is concluded that 
the concentration of both the pollutants increases as time 
increases and reaches their respective steady state values as 
t — > 00 . Due to the presence of sink phase the concentration of both 
the pollutants reduced considerably (compare Figures 2.4 and 2.5 
for 6j= 8 ^ = 0, 5j= §2 = 50 respectively). 

When flux is step function type at the source, the central 

line concentrations of both the pollutants decrease as time 

increases for t > t and the behavior is similar to the case of 

0 

constant flux for t ^ t^Csee Figures 2. 6-2. 7). Note that 
concentrations of both the pollutants have decreased significantly 
due to presence of sinks. 

The effect of variation of removal parameter associated with 
sink phase with vertical concentration profiles is shown in Figure 
2.8. It is observed from this figure that the concentrations of 
both primary and secondary pollutants decrease as removal 


parameter increases. 
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Finally, from the above discussion it may be concluded that 
the presence of sink phase in the second patch in the form of rain 
droplets, fog droplets or even externally introduced species 
(liquid or gas) can be very effective in removing the pollutants 
from the atmosphere. This analysis thus suggests a mechanism by 
which toxic gases leaked out in the atmosphere due to accidental 
discharge, etc. can be removed by introducing suitable liquid or 
gaseous phase in the environment using mechanical means. The above 
analysis can also be utilized for studies related to risk analysis 
particularly in cases of ^ accidental leakages of 
pollutants/toxicants in the atmosphere such as MIC leakage in a 


pesticide plant in Bhopal, India. 
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FIG.2.5 CONCENTRATION-DIST/MCE PROFILES FOR CONSTANT 
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FIG. 2. 8 VARIATION OF CONCENTRATION WITH REMOVAL 
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of dust released from an elevated source by taking into account 
the deposition and settling of dust particles under inversion 
condition. It may be noted that the sources are, in general, time 
dependent e.g., accidental explosion of a nuclear device, leakag^e 
of toxic material from storage tanks, etc. but very little 
attention has been paid to study the dispersion in such cases with 
gravitational settling, Cleary et a7.(l974), LI ewe 1 ynC 1983 ) . In 
particular, LI ewe 1 yn( 1983 ) presented the solution of unsteady 
state three dimensional transport diffusion equation with a first 
order reaction term for primary pollutant only emitted from a 
continuous point source using leakage at the base of stable 
inversion layer and deposition at the ground. 

Keeping in view the above, in this chapter we study the 
unsteady state dispersion of heavier air pollutant emitted from a 
time dependent point source forming secondary species in a two 
patch environment to see the effect of gravitational settling, 
which may be different in the two patches due to formation of fog, 
etc. The deposition at the ground and first order conversion of 
primary pollutant to secondary pollutant with removal of both the 
species has also been taken into account. To do so the unsteady 
state coupled three dimensional convective diffusion equations 
involving both primary and secondary pollutants have been solved 
incorporating the above mentioned factors. As in chapter II the 
different forms of time dependent flux at the source are 
considered in the analysis: (i) Instantaneous (ii) Constant, and 
Ciii) Step function type. 
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3 . 1 MATHEMATICAL MODEL 

Consider the unsteady state dispersion of heavier air 
pollutant, from a time dependent point source of strength Q(t) 
located at height h from the ground z = 0, forming secondary 

o 

species in a two patch environment under atmospheric inversion as 
shown in Figure 2.1 of the previous chapter. It is assumed that 
the time dependent source is located in the first patch and the 
wind is sufficiently large in both the patches so that downwind 
diffusion is negligible in comparison to advection. 

Under these assumptions the partial differential equations 
governing the concentration of pollutants in both the patches, 
with appropriate initial and boundary conditions can be written 
as follows, 


Patch I (0 ^ X < Xj , 0 i z H) 

The differential equation governing the concentration 
Cj(x,y,z,t) of the primary pollutant in the first patch is written 
as , 


dCj 

ar 


+ u 


ac 
1 aF 


1 


ac. 

^ = K , 
1 az yl 


a^c 


1 


a^c 


^ " ^zl 


1 


^ - (kj + kgj)Cj 


(3.1) 


ay" az 

where x, y, z are cartesian co-ordinates, t is the time, x-axis is 
taken in the downwind direction and z-axis vertically upwards, u^ 
is the mean wind velocity, w^ is the settling velocity which, 
however, depends on the size, shape and weight of the emitted 
particles but here it is assumed to be independent of position for 
simplicity, (Fisher and Macqueen, 1977). and are 
d i f f us i vi t ies , assumed to be constants, in y and z directions 


respectively. The constant k^ is the conversion rate coefficient 
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of primary pollutant to secondary pollutant and k., is its removal 
rate coefficient. 

The initial and boundary conditions for equation (3.1) are, 


^1 

(x ,y 

, z , t ) 

= 0 

at t 

= 0 

(3 

.2) 

^1 

(x ,y 

, z , t ) 

_ Q(t) 

■ 6(y) 

5(z-hg) at X = 0 

(3. 

3) 

^1 

(x ,y 

,z , t) 

= 0 

as y 

> ± 00 

(3. 

4) 


ac, 







V t 

zi az 

■ ^ ^ 

Cl =v 

c 

dl^l 

at z = 0 

(3 

.5) 


= 0 at z = H 


(3.6) 


where is the deposition velocity of the primary pollutant at 
the ground z = 0, H is the height of the inversion layer and 6( . ) 
is the Dirac delta function. 

Similarly the differential equation governing the 
concentration Cpj(x,y,z,t) of the secondary pollutant can be 
written as, 

ac 1 ac . . 

+ U - w = K — ^ + K — ^ + k.C - k .C . 
at 1 dx 1 az yl zl^^z 11 plpl 

(3.7) 

The initial and boundary conditions for equation (3.7) are, 
(assuming that there is no direct emission of secondary pollutant 
from the source) 


Cpi(x,y,z,t) = 0 at t = 0 


(3.8) 


Cpi(x,y,z,t) =0 at x = 0 ; y = 0 , z = h^ 


(3.9) 


Cpi(x,y,z,t) = 0 as y 


± 00 


(3.10) 
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K . + W,c 


z 1 dz 


. ^ = V , . C . at z 

1 pi dpi pi 


= 0 


(3.11) 


dC 


Ei=Oatz=H (3.12) 

is the deposition velocity of the secondary pollutant 


az 

where v 


dp 1 

at z = 0 and k , is its removal rate coefficient, 
pi 

Patch II (x 2: Xj, 0 ^ z ^ H) 

The differential equation governing the concentration 
C 2 (x,y,z,t) of the primary pollutant in the second patch is 
written as, 


ac, ac^ 
at^ ^2 ^ - 


dC^ 

2^ ^ ^ ^ 


d^C 

Z2 7^ - '“2 " « 

az 


where U 2 is the mean wind velocity, W 2 is the settling velocity, 
Ky 2 ^22 dif fusivities in y and z directions respectively. 
The constant k 2 is the conversion rate coefficient of primary 
pollutant to secondary pollutant and kg 2 is its removal rate 
coef f icient . 

The initial and boundary conditions for equation (3.13) are. 


C 2 (x,y,z,t) =0 at t = 0 


(3.14) 


C 2 (x,y,z,t) = Cj(x,y,z,t) at x = x^ 


(3.15) 


C 2 (x,y>z,t) =0 as y — > ± 00 


(3.16) 


ac 

^z2 dz~ '^2S ~ '^d2‘^2 


at z = 0 


(3.17) 


ac, 

d 

W 


= 0 at z = H 


(3.18) 


where v ^2 is the deposition velocity of the primary pollutant at 
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z = 0 . 


Similarly 
concentration 
as , 


the 

Cp2(x.y.z 


differential equation governing the 
, t) of the secondary pollutant is written 


ac 


P2 


ac 


at 


+ u 


P2 


2 ax 


w - V 

2 az y2 


a^c 


P2 ^ 


a^c 


K 


P2 


ay 


az^ 


V 2 ~ ^p2^p2 


(3.19) 


The initial and boundary conditions for equation (3.19) are, 


C „(x,y,z,t) =0 at t = 0 

p2 


Cp 2 (x,y,z,t) =Cpj(x,y,z,t) at x=Xj 


Cp2(x,y,z,t) =0 as y — » ± 00 

ac 2 

^z2 ai^ ^ "'2S2 = '^dp2S2 at 2 = 0 


ac 


p2 


az 


0 at z = H 


(3.20) 

(3.21) 

(3.22) 

(3.23) 

(3.24) 


where v^p 2 is the deposition velocity of the secondary pollutant 
at z = 0 and k ^ is its removal rate coefficient. 

p2 

As in the previous chapter, the following forms of flux Q(t) 
are considered in the subsequent analysis, 

(i) Flux is Instantaneous 

Q(t) = 6(t) 

(ii) Flux is Constant (3.25) 

Q(t) = 

(iii) Flux is Step function type 

Q(t) = Q^, 0 i t :S t^ 
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To cast the problem in dimensionless form, the following non- 
dimensional quantities are defined, 


K t 

t = — 12_, X 


u 


K X 
zo 


u H 
0 


2 ’ 


y = 


H ’ 


z = 


H 


^’i u 


u H 

C. = -7^ C. , 
1 Q 1 


u H 

C = C 

pi Q. pi 


h = -tr 
s H 


QCt) = 


Q, 


- _ Qo ^zo 


H w 


w . 
1 


K 


zi 


1.2 


where K and u are reference diffusion coefficient and wind 
zo 0 

velocity respectively. 

The equations (3.1-3.25) can be written in dimensionless 
form (dropping bars for convenience) as: 


Patch I (0 5 X < x^, 0 ^ z ^ 1) 

ac. ac. ac. a^c. a^c. 

_ + Uj - Wj ^ ^ + 7^ — -J- 

ay az 

Cj(x,y,z,t) =0 at t = 0 
C.(x,y,z,t) = 6(y) 6(z-h ) at x 

X ^ ^ 

Cj(x,y,z,t) =0 as y — > ± oo 

ac. 

az~ '^l^l ~ ^1^1 at z = 0 


(a, 




0 


(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 


az 

ac 


= 0 


pi 


at 

ac 


z = 1 


(3.31) 


at 


+ u 


1 ax 


pi _ 


ac 


w 


pi _ 


a^c 


1 az 


= ^ 


' ay2 


Pi ^ 


a^c j 

rT— * “1=1 

dz 


“pi Si 


(3.32) 
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Cpi (X.y ,z, t) 

= 0 

at 

t = 0 


(3.33) 

Cpi(x,y,z,t) 

= 0 

at 

X = 0; y = 0 , 

z = h 

s 

(3.34) 

Cp j ( X , y , z , t ) 

= 0 

as 

y — > ± 00 


(3.35) 

ac , 
ai^ * 

= N .C , 
pi pi 

at z = 0 


(3.36) 

ac , 

az” = ° 

Z 

= 1 



(3.37) 

Patch II (x 

^ ^1 

, 0 

^ Z 1') 



ac^ ac2 

at ^^2 a5r ■ 

^2 

ac2 

az 

a^c. 

= ^2 — 2 ^ ^ ^2 
^ ay^ ^ 

a^C2 

- '“2 - °g 2 )S 

(3.38) 

C2Cx,y,z,t) = 

0 

at 

t = 0 


(3.39) 

C2(x,y,z,t) = 

Cj Cx ,y , 

z , t ) at X = 


(3.40) 

C2(x,y,z,t) = 

0 

as 

y — > ± 00 


(3.41) 


dC^ 

^ ^ ^2^2 = ^202' at z = 0 (3.42) 

dC^ 

= 0 at z = 1 (3.43) 


dC 


p2 


ac 


at 


+ u 


P2 


2 ax 


Cp2(x,y,z,t) = 
Cp2(x,y,z,t) = 


ac 


— w 


2 az 


p2 _ p 

P2 


a^c 


P2 + 


a^c 


P2 


ay 


'2 az2 


“ 2 S 


“P 2 S 2 


O at t = 0 


(3.44) 

(3.45) 


C .(x,y,z,t) at X = X. (3.46) 

pi 1 


Cp2(x,y,z,t) 


O as y — > ± CO 


(3.47) 
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ac 

az '^2^p2 " ^p2‘^p2 z = 0 


ac 


az 


= 0 at z = 1 


C3.48) 


(3.49) 


k.H‘ 

where a. = 

1 K 

zo 
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gi 


tgiH 

■‘zo 


K . 
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zi 
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z 1 

K 
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pi 


V , . H 

dpi 


z 1 


“pi 


, (i =1,2) 


k 

pi 

~K 

Zo 


(i) Q(t) 

(ii) Q(t) 

(iii) Q(t) 


Q 5(t) 
0 


= 1 
= 1 , 
= 0 . 


0 ;S t i t 


t > t 


(3.50) 


3.2 METHOD OF SOLUTION 

Equations (3.26) and (3.32) can be written in compact form 

as , 


fc:, n 


1 


c . 


L pij 



“l “gl 


~ a. 


0 


a 


pi 


'1 


^Pl 


= 0 


(3.51) 


where is the advection/dif fusion operator 


Z s - — 

^1 at 


^1 fe ^ ^1 ii ^ ^1 


a^ ^ a2 

— •" ^1 — 


ay az“ 

Following the same method as in the previous chapter, we get 
the uncoupled system for patch I, 


- (ttj + ~ ® 


(3.52) 


^1®1 " “pl®l = 


0 


(3.53) 
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“l 

where B, = C , - — ^ C, (3.54) 

1 pi «p, - a, - 1 

The conditions for remains same. The corresponding 

conditions for are, 

Bj(x,y,z,t) =0 at t = 0 (3.55) 

B (x,y,z,t) = - 5(y) 5(z-h ) at x = 0 (3.56) 

pi “l “gl 1 ® 


Bj(x,y,z,t) =0 as y — > ± oo 


(3.57) 


aB ' N - N, 

^ 

dz pi 11a 


‘pi “l 


a, C. 

“gi ‘ ' 


at z = 0 


(3.58) 


aB 


1 


az 


= 0 


at z = 1 


(5.59) 


Similarly we also get the uncoupled system corresponding to 
equations (3.38) and (3.44) for Patch-II, 


^ 2*^2 ^“2 “ g 2^^2 “ ® 


^2®2 - «p2 ®2 = ° 


(3.60) 

(3.61) 


a 

m: 


a a 

‘2 ax '^2 az 


£ 


where e - ^ - u^ ^ ^ ^2 

dz 


(3.62) 


a. 


and = C — 

2 p2 a 


p2 - “2 ■ “g2 2 


(3.63) 


The conditions for C 2 remains same. The corresponding 
conditions for B 2 are, 


B2(x,y,z,t) = 0 


at t = 0 
r 


(3.64) 


B 2 (x,y,z,t) = Bj(x,y,z,t) + 


a, 


a. 


“pr“r“gl “p2"“2"“g2 j 

at X = X, 


(3.65) 
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B 2 (x,y,z,t )=0 as y — >±® 


(3.66) 


6 B, 

i . 

w 


(N „ - w„)B., + — 
p 2 2 2 a 


- N, 

P 2 A 


p 2 “2 


a 


g 2 
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z = 0 


(3.67) 


aB, 

i 

W 


= 0 at z = 1 


(3.68) 


Introducing the following transformations 


w, 

Cj(x,y,z,t) = C^^(x,y,z,t) exp j - 


w 

(z-h ) - — 

s 47 


— ’‘I 


(3.69) 


Bj^(x,y,z,t) = Bj^(x,y,z,t) exp 
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w . 


- 75 — (z-h ) 
2 Yj s 


4 YjUj 


X (3.70) 
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4 J 


(3.71) 


B 2 (x,y,z,t) = B 22 (x,y, 2 ,t) exp 


w. 


2 Y. 


(z-h ) - 


^2 'I 

4 


(3.72) 


Equations (3.52 - 3.68) can now be written as, 


Patch I (0 ^ X <Xj , 0 £ z ^ 1) 
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Patch II(xkXj, 0;SZ^1) 
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3.3 CONCENTRATION DISTRIBUTIONS 

Solutions for and ^22’ ®22 obtained using 

Laplace and Fourier transform techniques, (Carslaw and Jaeger, 
1959) in each cases as, 

3.3.1 CASE I WHEN FLUX IS INSTANTANEOUS AT THE SOURCE 
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where 5(t - x) is the Dirac delta function defined in chapter II. 
And , 
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^2 

G (z) = cos u z + — sin u z 
n 


N 

H^(z) = cos ^ z + sin F z 

n ^n F ^n 

n 


N 

J (z) = cos p z + -i— sin p z 
n n p n 

n 


N. 


r 1 7 

P = (cos X„z + T — sin X z) dz 
n J n X n 

0 ” 

^2 2 
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n J ^n u ^n 
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P_„ = J (cos p_z + 77 ^ sin p_z)'“dz 


n2 
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_I 

n p. 


0 


n 


n 


^n’ ^n’ ^n roots of the following transcendental 

equations respectively, 
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3.3.2 CASE II WHEN FLUX IS CONSTANT AT THE SOURCE Q(t) = 1 

C3.101) 


'11 


(x,y,z,t) = — ^ PCx,y,z) H^^t - 


Bjj(x,y,z,t) = Q(x,y,z) H[^t - 


(3. 102) 


CjjCx.y.z.t) = -i- M(x,y,z) h[i " (-nj * "S^]] 


(3.103) 


B 


22 


(x,y,z,t) = N(x,y,z) - (-ji + 


(3. 104) 


From equations (3 . 101 )-{3 . 104) it is observed that if t — > oo 
the expressions reduce to the steady state solutions. 

3.3.3 CASE III WHEN FLUX IS STEP , FUNCTION TYPE AT THE SOURCE 
Q(t) = 1 , 0 i t ^ t^ 

= 0 , t > t 


Cjj(x,y,z,t) = 4- P(x.y,z)[H(t - (3. 105) 
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j (x ,y ,z , t) = 


4- C!(x,y,z)[H[t - ^]H(t-t^)](3.106) 


y.z.t) = 4- H(x,y,z) [H[t -(4i * 

- - ‘o - 


C22'’‘- 


(3.107) 


B 22 (x,y,z,t) = -L. N(x,y.z) [hJi 

't - t - )] 

0 I Xi^ )) 0 j 


(3.108) 


where H(t-x) is the Heaviside function defined in chapter II. 

From equations (3 . 105)-(3 . 108) it is seen that as t — ¥ oo, 

0 

these equations reduce to equations ( 3 . 1 01 )- ( 3 . 1 04 ) 

respectively. 

The final solutions for , C 2 . B 2 can be obtained from 

equations (3.69)- (3.72) by putting the values of Cjj. ®1 1 ’ *^22’ 
®22 above equations in each cases. The dimensionless 

concentration and 0^2 can then be obtained from , B^ and C 2 , 

B 2 in each of the above cases by the following relations, 


Cpi(x,y,z,t) = Bj(x,y,z,t) + — 


pi 


“1 “gl 


OCr 


C ^(x.y.z.t) = B„(x,y,z,t) + 

p2 2 ap2 - «2 " “ 


g2 


Cj (x , y , z , t ) 


C2(x,y,z,t) 


3.4 DISCUSSION AND RESULTS 

To see the effect of various parameters the expressions for 
concentration distribution of both primary and secondary species 
in both the patches have been computed and displayed graphically 
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in Figures 3.1-3. 12. Tlie gjaa'amie ters ^sal'Ciiies used in the analysis 
are aj-a2-2.77, “«g2 " ^ ^ 20.0, Nj = 4.0, N2=4.2, 

Npi=0.40, Np2=0-42, vVj&w^ = 0 .0 ,1,5, ■uiii|=*W2 = 1.0, Pj=p2 “ 10.0, 

y =)- = 1,0, h =0.20, x,=0,03 , y=0-0 , Q 0 , t =0.04, e=0,005. 

It is observed that the ccnnceMnt iiions of both primary and 
secondary pollutants decxe as e a:s doownwroniid distance increases in 
all the three cases axd as lth»,« rsat e ‘Cooff conversion of primary 
pollutant to secondary p>ol lLjtaa“t imt riessiisj^jes the concentration of 
primary pollutant decreas ess amd thuatt of secondary pollutant 

increases at a particular t ime and 1 ocaiattii 0);n. 

When flux is i ns tantame-ou s ad; tbes aeource, the concentration 
distributions of both "the pol lu ta_nts Wt'Wh downwind distance are 
shown in Figures 3. 1-3. 2 fcor d:if:ie resit *<*al!ues of time t along z = 
0.05 (below the source h«i gh_t> a, niffi! & W 2 = 0.0, 1.5 

respectively. It is seen tiha t tine csioncs'Wit rations of both primary 
and secondary pollutants decr'essoe assttiDes increases. On comparing 
these figures it is also n«t «d that f!<o»in' the points located belo-w 
the source height, the ooncenttrat iocis (^IfK pollutants increase as 
the settling velocity increas ess, How«ev'err',, for the points above the 
source height, the con cent rst:io i de acp'Ssases as the settling 
velocity increases, (See ILgiures 3 .8J i^Uwere we have shown the 

vertical concentrati on prof il esi cof pirl pollutant for different 

settling velocities in the <a se oof i«ns IspBvrtltaneous flux). 

When flux is constant ait tl he sso urcioesj , the concentrations of 
both the pollutants with downwiincsi dL-5’t)i:iD'ffli!t8 are plotted in Figures 
3.3 and 3.4 for different vaEl”ue=s of tti nei ttt along z = 0.05 and & 
W 2 = 0.0, 1.5 and it is se*ea Bhaat tttu'e concentrations of these 
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pollutants increase as time increases and reaches their respective 
steady values as t — > 00 . On comparing the figures it is again 
noted that for the points located below the source height, an 
increase in settling velocity increases the concentrations of 
pollutants and reverse is the case for the points above the source 
he i ght . 

When flux is step function type at the source, the 

concentrations of both the species decrease as time increases for 

t > t and behavior is similar to the case of constant flux for t 
0 

t . The increase in settling velocity increases the 
0 

concentrations of both the species. (See Figures 3.5 & 3.6 ). 

To see the effect of settling velocity more clearly we have 

also shown horizontal and vertical concentration profiles at a 

fixed time for constant and instantaneous sources. 

In Figure 3.7 the concentration profiles of primary pollutant 

with downwind distance are shown for different values of settling 

velocities w^ and ^ fixed time t = 0.06 at z = O.OCthe 

ground) and z = 0.5 (above the source height) for constant flux 

at the source. It is quite apparent from the figure that with the 

increase in settling velocity, the concentration of the pollutant 

increases at the ground level while it decreases at z = 0.5 i.e., 

for the points above the source height h^ = 0.2. 

Similarly in Figures 3.8 and 3.9 we have shown vertical 

concentration profiles for different settling velocities at 

different distances x for instantaneous and constant flux 

respectively. Note the concentration increase with increase in 

settling velocity below the source height h = 

s 


0.2 


and 
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concentration decrease above the source height. 

In Figure 3.10 the vertical concentration profiles of 
secondary pollutant are drawn for constant flux at the source and 
similar observations are made. 

From the above analysis one important thing has been observed 
that as settling velocity increases, the concentration of 
pollutants above the source height decreases but increases below 
the source height. This is expected because after emission of 
pollutants from the source gravity comes into effect and heavier 
particles are likely to have higher settling velocity in 
comparison to not so heavier particles. 

The effect of variation of ground level deposition velocity 
on the concentration distribution of primary pollutant with 
vertical height for a constant source at different distances x are 
shown in Figures 3.11 and 3.12. It is seen that the concentration 
decreases with increase in deposition velocity and this effect is 
more predominant below the source height. Also on comparing 
figures 3.11 (zero settling velocity, '^ 1=^2 = 0.0 ) and 3.12 (w^ 
and w„ = 1.0 ) it is seen that the concentration has increased 
considerably below z = 0.2 due to increase of settling velocity 
but has decreased above z = 0.2. 



DIMENSIONLESS CONCENTRATION 



DIMENSIONLESS DOWNWIND DISTANCE 

FIG. 3,1 CONCENTRATION-DISTANCE PROFILES ALONG 2=0.05, 
WITHOUT SETTLING ( INSTANTANEOUS FLUX ). 
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DIMENSIONLESS DOWNWIND DISTANCE 


FIG.3.2 CONCENTRATION-DISTANCE PROFILES ALONG z=0.05, 
WITH SETTLING ( INSTANTANEOUS FLUX)- 
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DIMENSIONLESS DOWNWIND DISTANCE 

FIG.3.3 CONCENTRATION-DISTANCE PROFILES ALONG z=0.05 
WITHOUT SETTLING (CONSTANT FLUX)- 




DIMENSIONLESS CONCENTRATION 


90 


0.4 


— Primary 

... secondory P®"" 

, = 0.05 

V =0-0 

^ .1.5 
!!'' -5 
r.^;o.o3 


■♦= 0.02 ,, 0 ^^ 1=006 


Q03 


0.05 


0.07 


Kin distance 
DIMENSIONLESS DOWNWIN 

IG. 3. 4 CONCENTRATION-DISTANCE ). 

WITH SETTLING (CONSTAN* 


0.05 



- Primary pollutant 
" Secondary pollutant 


91 



UJ 

u 

2 

< 


Q 

Q 

2 


2 


O 

Q 

cO 

to 

LU 


o 

to 

LU 


2 

Q 


N0liVdiN3DN0D SS31N0ISN3Nia 


FIG.3.5 CONCENTRATION- DISTANCE PROFILES ALONG z = 0 05 WITHOUT 
SETTLING (STEP FUNCTION TYPE FLUX ). 
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DIMENSIONLESS CONCENTRATION 
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DIMENSIONLESS DOWNWIND DISTANCE 

FIG,3.7 EFFECT OF SETTLING VELOCITY ON CONCENTRATION 
DISTRIBUTION FOR CONSTANT FLUX. 
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i VERTICAL CONCENTRATION PROFILES FOR DIFFEf^h 
SETTLING VELOCITIES ( INSTANTANEOUS FLUX ). 




DIMENSIONLESS VERTICAL HEIGHT 


FIG. 3. 
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DIMENSIONLESS CONCENTRATION 

VERTICAL CONCENTRATION PROFILES FOR DIFFERENT 
SETTLING VELOCITIES ( CONSTANT FLUX). 



DIMENSIONLESS VERTICAL HEIGHT 



FIG.3.10 VERTICAL CONCENTRATION PROFILES OF SECONDARY 
POLLUTANT FOR DIFFERENT SETTLING VELOCITIES 
( CONSTANT FLUX ). 




DIMENSIONLESS VERTICAL HEIGHT 



DIMENSIONLESS CONCENTRATION 

FIG.3.11 VERTICAL CONCENTRATION PROFILES FOR DIFFERENT 
DEPOSITION VELOCITIES, WITHOUT SETTLING (CONSTANT 
FLUX ). 
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FIG.3.12 VERTICAL CONCENTRATION PROFILES FOR DIFFERENl 
DEPOSITION VELOCITIES; WITH SETTLING ( CONSTANT 
FLUX). 




CHAPTER IV 


Unsteady State Dispersion of a Reactive Air Pollutant from a Time 
Dependent Point Source : Effect of Greenbelt Plantation 

4 . 0 INTRODUCTION 

The pattern of pollutant dispersion is modified by the 
various removal mechanisms. One such mechanism is the plantation 
of greenbelt, between the source and the receptor, (e.g., a 
habitat, a monument) capable of absorbing pollutants and inducing 
deposition at the surface of the canopy of trees. It is well 
established that certain plants have the capability of removing 
pollutants by absorption, deposition, etc. It may be mentioned 
here that these characteristics of greenbelt plantation not only 
depend upon the distribution, size, height and nature of the 
canopy but also on the type of pollutant, .species of trees, 
shrubs, herbs, etc. planted in the greenbelt. Thus, the choice of 
the plantation should be such that the species planted be 
relatively resistant to the exposure of the specific pollutant. 
Number of investigations have been made in the past to study the 
effect of pollution on plant species, e.g. DainesC 1968) , 
TreshowC 1968) , Davis(1972), Takemoto and Noble(1982), 

Dochinger ( 1988) , Evans(1988). In particular, survey of Indian 

plants in relation to atmospheric pollutants have been conducted 
by Ahmad et af.(1988) who have recommended a number of pollutant 
tolerant trees, shrubs and herbs for landscaping urban and 
industrial areas for mitigating air pollutants. The modification 
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of dispersion pattern of gaseous and particulate matters within 
plant and vegetative canopies has also been studied by many 
authors, e.g., MacBean( 1 968 ) , Hill(1971), Bennet and Hill(1973, 
1975), ChamberlainC 1975a, b) , O’Dell et 37.(1977), Petit et 
a/. (1976), Bache(1979a,b) , Smith(1981), Hosker and Lindberg ( 1 982 ) , 
SI inn( 1982a) , Kapoor and Gupta(1984), Wiman and Agren(1985), Viman 

and Lannef ors( 1985) , Wiman et ai. (1985), Baldocchi et a7.(1987), 
Shukla and Chauhan( 1 987 ) , Shukla et a7.(1987). In particular, 
Petit et a7.(1976) presented some results concerning 
characteristics of air flow within and above a forest and 
calculated SO 2 fluxes at the top of the canopy to understand the 
diffusion mechanism of SO 2 . Bache(1979) has also developed some 
analysis and suggested a modified form of diffusion equation to 
understand particulate transport mechanism within plant canopies 
and also compared analytical solutions with numerical results. 
SI inn( 1982a) gave a theoretical framework to predict particle 
deposition to vegetative canopies and obtained approximate 
so tions of diffusion equation describing the air concentration 
of particles using variable wind velocity profile. A review of the 
atmospheric deposition and plant assimilation of gases and 
particles has been presented by Hosker and Lindberg( 1982 ) and 
Smith(1981). Wiman and Agren(1985) developed a two dimensional 
steady state model describing aerosol concentration within forest, 
assumed to be very long and homogeneous in the direction 
perpendicular to the wind direction, by considering the 
interaction between forest structure, open field, forest 
aerodynamics and aerosol characteristics. 
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From these studies it may be speculated that if a greenbelt, 
capable of absorbing the pollutant, is planted suitably between 
the source and the region to be protected, it is possible to 
minimize the effect of pollution on this region. F 1 emm ingC 1 967 ) , 
BernatzkyC 1968) , Hanson and Throne ( 1 970) , WarrenC 1973) , Gupta and 
Kapoor ( 1985) . In this direction Kapoor and Gupta(1984) have 
studied the effect of greenbelt on reduction of concentration of 
the pollutant emitted from ground level sources under steady state 
conditions using Gaussian plume type model. 

It may be noted here that in real situations, however, 


pol lutant 

sources are time dependent (such as 

in the case 

of 

accidental 

explosion of toxic 

gas storages) , it 

is necessary 

to 

model the 

phenomenon, in 

general , as 

a dynamic process which 

includes 

the effects of 

stack 

height , 

1 ocat ion , 

thickness 

and 

height of 

the greenbelt, 

the 

location 

and the 

height of 

the 

receptor , 

etc . 







In this chapter we, therefore, study the unsteady state 
dispersion of a reactive air pollutant emitted from a time 
dependent point source and its reduction due to greenbelt. The 
optimum location of the greenbelt between the source and the 
receptor for maximum concentration reduction of the pollutant at 
the receptor has also been studied. The problem is investigated by 
solving the unsteady state three dimensional convective diffusion 
equation with removal terms using suitable matching and boundary 
conditions . 
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4.1 MATHEMATICAL MODEL - I : THE CASE WITH DEPOSITION AT THE TOP 
OF THE GREENBELT 

Consider the unsteady state dispersion of a reactive air 
pollutant from a time dependent point source in presence of 
greenbelt under inversion condition. It is assumed that the 
greenbelt of mean thickness d and height h is located at a 

5 

distance x, from the source of height h in the wind direction 
extending in y direction symmetrically . The physical situation is 
illustrated in Figure 4.1. In modelling the process of the effect 
of greenbelt on pollutant distribution, we are ignoring the 
ecological aspects such as the internal mechanisms like 
distribution of leaves, locations of deposition, effect of 
pollutant on plants, etc. as our main concern here is to see the 
effect of gross depletion parameter characterizing the greenbelt 
which is acting as a sink. Petit et ai. (1976), Bache(1979), 
SI inn( 1982a) . In the greenbelt region we use simplified convective 
diffusion equation with removal term characterizing the greenbelt 
in such a way that different removal parameters can be chosen for 
different greenbelts, SI inn( 1982a) , Wiman and AgrenC1985). The 
vertical and horizontal diffusion processes in this region are 
also neglected due to denseness of the plant canopies, shrubs 
grown on the ground below canopies, etc. However, in the regions 
without greenbelt, usual unsteady state convective diffusion 
equations with reaction terms are considered with suitable 
matching and boundary conditions at the interfaces surrounding the 
greenbelt. In particular, the condition used at the top of the 
greenbelt (i.e., at the interface of the greenbelt region and the 
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region above it, the region III) implies that the deposition flux 
is proportional to the difference between the concentration of the 
same pollutant within and outside the greenbelt. Region IV is the 
region where receptor is located. It is also assumed that the wind 
is sufficiently large in the regions I, III and IV so that the 
downwind diffusion is negligible in comparison to advection. 

Under these assumptions the dynamic equations governing the 
concentration distribution of the pollutant in all the regions can 
be written as, 


Region I (0ix<Xj , 0^Z:sH) 

The differential equation governing the concentration 


Cj(x,y,z,t) of the pollutant in the region I is written as. 


ac^ 

at 


ac 


+ u 


1 


a^c, 


1 ax 


= K 


ay 2 


a^c 

" Si - ss 

az 


(4.1) 


where x, y, z are cartesian co-ordinates, t is the time, x-axis is 

taken in the downwind direction and z-axis vertically upwards, u^ 

is the mean wind velocity, K , and K , are dif fusivities in y and 

yl zl 

z directions respectively, taken as constants. The constant k^ is 
the reaction rate coefficient of the pollutant. 

The initial and boundary conditions for equation (4.1) are. 


C^(x,y,z,t) =0 at t = 0 


C.(x.y,z,t) = 5(y) 5Cz-h^) 

X ^ 


CjCx,y,z,t) =0 as y — * ± a> 


at X = 0 


(4.2) 

(4.3) 

(4.4) 


ac. 

K , ^ = V..C, 

zl az dl 1 


at z = 0 


(4.5) 
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aCj 

0z~ 


= 0 


at z = H 


where is the deposition velocity of the pollutant at z 

5( . ) is the Dirac delta function. 


(4.6) 
0 and 


Region II (x^ ^ x< 0 i z < h^) 

The differential equation governing the concentration 
C 2 (x,y,z,t) of the pollutant in the region II is written as, 


ac, 

dt 


•4" U 


dCr 

d 

2 W' 


a^c. 


= K. 


ay2 


(k2 + X)C2 


(4.7) 


where U 2 is the mean wind velocity inside the greenbelt, Ky .2 
K ^2 d i f f us ivi t i es in y and z directions respectively, k 2 is 

the reaction rate coefficient of pollutant and X is the depletion 
rate coefficient due to greenbelt which is acting as a sink. 

The initial and boundary conditions for equation (4.7) are. 


C 2 (x,y,z,t) =0 at t = 0 


(4.8) 


C 2 (x,y,z,t) = Cj(x,y,z,t) at x = 


(4.9) 


C 2 (x,y,z,t )=0 as y 


± <B 


(4.10) 


Region III (x^, ^ x < X 2 > h^ i z H) 

The differential equation governing the concentration 
C 2 (x,y,z,t) of the pollutant in the region III is written as, 


^=3 ^ „ ^=3 . „ 

at ^3 ax y3 _ 

ay 


a^c. 

^ " ^Z3 ^ 
0Z 


- ^3^ 


(4.11) 


where u^ is the mean wind velocity, and dif fusivities 
in y and z directions respectively. The constant k^ is the 
reaction rate coefficient of the pollutant. 

The initial and boundary conditions for equation (4.11) are, 
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C 2 (x,y,z,t) = 0 

at t = 0 

(4.12) 

Cglx ,y , z , t ) = 

(x,y,z,t) at X = Xj 

(4.13) 

C 2 (x,y,z,t) = 0 

as y — ^ ± OD 

(4 .14) 

ac„ 

^z3 az ''d3^‘^3 

- C„) at z = h . 

2 g 

(4.15) 

ac.. 

az = ° ^ 

= H 

(4 . 16) 

where v^^ is the deposition velocity 

of the pollutant at z = h . . 

Here equation(4. 

15) implies that due 

to presence of greenbelt, the 


deposition is taking place at the top of the greenbelt (z=hg) and 
the corresponding flux is proportional to the difference of 
concentration of the same pollutant within and outside the 
greenbe 1 t . 


Region IV ( ( x i X 2 5 0 :S z H) 

The differential equation governing the concentration 
C^(x,y,z,t) of the pollutant in the region IV is written as, 


dc ac a^c a^c 

^ " "4 -^ = ^y4 -r/ " ^Z4 -7:^ - ^4^4 


(4. 17) 


ay' 


az' 


where u^ is the mean wind velocity, and are dif fusivities 
in y and z directions respectively. The constant is the 
reaction rate coefficient of the pollutant. 

The initial and boundary conditions for equation (4.17) are, 


C^(x,y,z,t) =0 at t = 0 (4.18) 

C^(x,y,z,t) = C 2 (x,y,z,t) at x = X 2 ; z < h^ 

= CgtXjy.z.t) at X = X 2 ; z ^ h^ 


(4.19) 



106 


C^(x,y,z,t)=0 as y — >±cd 


ac 

I . - — = Vj-C. at z = 0 
z4 dz d4 4 


(4.20) 

(4.21) 


dC. 

^ =0 at z = H (4.22) 

dz 

where is the deposition velocity of the pollutant at z = 0. 

As in the previous chapter the following forms of time 
dependent flux Q(t) are considered in the analysis: 

(i) Flux is Instantaneous, Q(t) = 5(t) 

(ii) Flux is Constant, Q(t) = Q (4.23) 

o 

(iii) Flux is Step function type 

Q(t) = Q , 0 t ^ t 

C 0 

=0, t > t 

0 

Define the following dimensionless quantities, 


K t 

T = 


^zo^ - y - 2 - ^s - ^i 

X = y = -ff-" z ^s ’ ^i ^ iT 


(i = 1,2) 


^i = 


u H 
0 

«. 


Cj, Q(t) = 2 ^ 


7=f _ Qo 20 


X = 


H^X 


"C c 

where K and u are reference diffusion coefficient and wind 


zo 0 

velocity respectively. 

The equations (4,l)-(4.23) can be written in dimensionless 
form (dropping bars for convenience) as : 


Region I (0ix<x., Oizil) 


ac ac. 

at- '^i -w = ^1 


a^c, 


a^c. 


ay 


2“ ^1 


dz 


ir “i^i 


(4.24) 
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Cj(x,y,z,t)=0 at t=0 


C, (x.y.z.t) = 6 (y) 5(z-hg) at x = 0 


Cj(x,y,z,t) =0 as y — » ± co 


ac 

— - = N.C. at z = 0 
az 11 


ac j 


= 0 at z = 1 


Region II (x,i x <x,. 0 z <hg) 


ac. 


ac. 
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^ ^2 ax ■ ^ 
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2 ay2 2 2 


C 2 (x,y,z,t) =0 at t = 0 


C 2 (x,y,z,t) = Cj(x,y,z,t) at x = x^ 


C 2 (x,y,z,t) =0 as y — > ± oo 
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CgCx.y.z.t) =0 at t = 0 


CgCx.y.z.t) = Cj(x,y,z,t) at x = x^ 
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± 0t> 


ac, 
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= 0 at z = 1 


= h 
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zo 


N. = 
1 
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(4.45) 
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4.2 SOLUTION OF THE MODEL 

The solution of equations (4.24-4.45) are obtained using 
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in each cases as follows, 
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4.2.1 CASE I WHEN FLUX IS INSTANTANEOUS AT THE SOURCE 
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4.2.2 CASE II WHEN FLUX IS CONSTANT AT THE SOURCE Q(t) = 1 
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4.2.3 CASE III WHEN FLUX IS STEP FUNCTION TYPE AT THE SOURCE 
Q(t) = 1 , 0 t i 

=0. t 
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4.2.4 DISCUSSION 

To discuss the effect of various parameters on the 
concentration distributions of pollutant, the expressions obtained 
above have been computed, choosing suitable dimensionless values 
and displayed graphically in Figures 4.2 - 4.17. The parameters 
values chosen are Q = 1.0, u. = 1.0, (i = l,3,4), u., = 0.05, y. = 

0 1 Z 1 

1.0, (i=l,3.4), a. = 2.77, \ = 10.0, = 0.5, = 4.0, = 0.4, 

p. = 10.0, Ci = l,3,4), = 0.5, location of receptor x = 0.1, y = 

1 

0.0, width of greenbelt d = 0.002, height of the greenbelt h = 

o 

0.05, t = 0. 12, e = 0.005. 

0 

It may be noted here that the height of the source plays an 

important role in determining the optimum location of greenbelt 

between the source and the receptor to affect the maximum 

concentration reduction. We have, therefore, considered three 

cases viz, h > h , h = h , and h < h . 

s g s g s g 

When flux is instantaneous at the source, the concentration 
distributions with and without greenbelt at x = 0 . 1 (receptor) in 



chansing the greenhe.t location x, 
,eg.on IV are P eoncentrat ion reduotton 

towards the receptor ,2-44 respectively at 

tor h, = 0.., O.OS, 0.0. in «.vres 4.^ ^ 

different receptor "^'^ 2 ^''.. o,05) and at time t = 0.102. 
above the greenbelt hei g t is seen that maximum 

When h^ > hg CFigure. 4.2) ' greenbelt 

concentration reduction takes concentration reduction 

.s .coated cioser to the receptor, and h 
is greater if the height of the recept 

, .. i, (Figure 4.3) and h, < (Figure 

For the case when , g 

j +v,a+ the greenbelt snouts 

,. 4 ,, it is tcnnd that 

source to affect maximum concen 3 „„rces are ground 

Tt is therefore, speculated that if the 

, c, it is desirable to surround this area by 
based lihe storage tanhs, ,„iioiants. It is also 

>,»lt to avoid further spreading of P 
the greenbelt ooncentration reduces 

noted from these figures that 

considerably due to greenbelt. concentration 

Vhen flux IS oonstan receptor x = 0.1 

.iatrihutions„ithand»ithout.eenbe.ta^^^^^^^^ 

are plotted m Figures heights 

location x, towards the receptor ^ ^ 

as in the ^ this case also similar results are 

heights z = 0.02 and 0.09. In 

noticed as before. source, the 

is step function type at 

When flux IS acres 4.8 - 4.10. It is 

f ie-c are plotted m Figur 

concentration profiles 


noted that for t i t^ the behavior is same as in the case of 

constant flux but when t > t the concentration becomes zero at 

0 

the receptor x = 0.1. 

In Figures 4.11 - 4.13 vertical concentration profiles at the 

receptor location x = 0.1 are shown for different locations of 

greenbelt for three cases of source heights as before for constant 

flux only. It is again found from these figures that for the 

maximum concentration reduction for h > h_ , the greenbelt should 

s g 

be nearer to the receptor while for h^ i h^ , the greenbelt should 

be closer to the source as pointed out earlier. 

In Figure 4.14 the concentration profiles with and without 

greenbelt with downwind distance are shown for h = 0.2, h^ = 

s g 

0.05, d = 0.004, t = 0.2 at two values of z = 0.02, 0.09. It is 

seen that the concentration is same with and without greenbelt 
before the greenbelt but in the greenbelt region x^ = 0.06, d = 
0.004, concentration decreases significantly and then increases 
slightly beyond the greenbelt and there after decreases 
continuously for the points away from the greenbelt. It is 
concluded here that if the receptor is far away from the greenbelt 
the concentration difference between with and without greenbelt 
approaches zero. 

In Figure 4.15 the effect of depletion coefficient X due to 

greenbelt is shown on the concentration distribution at the 

receptor x = 0.1 for h = 0.2 for constant flux and it is 

s 

found that the concentration decreases as X increases. 

In Figure 4.16 the concentration profiles at the receptor x = 
0.1 for constant source are shown for different deposition) 



velocities at the top of the greenbelt located at = 0.08 for h^ 
= 0.2,. It is seen that as the deposition velocity increases the 
concentration at the receptor decreases. 

In Figure 4.17 the vertical concentration profiles at 
different receptor locations x = 0.06, 0.10, 0.14 are displayed 
while the greenbelt is located at x^, = 0.04. It is found that as 
the receptor moves away from the greenbelt, the concentration 
decreases at the receptor. Also the concentration difference 
between the cases with and without greenbelt, approaches zero if 
the receptor is far away from the greenbelt. 


4.3 MATHEMATICAL MODEL - II : THE CASE WITH INTER - EXCHANGE 

PROCESS IN THE ADJACENT LAYER TO THE TOP OF THE GREENBELT 

In this case we again consider the atmosphere to be divided 

into four zones between the source and the receptor as in case 

4.1. We consider here the inter-exchange process in the layer of 

thickness H adjacent to the top of the greenbelt in the region 

III. We assume that in this layer the decrease in the rate of 

aCg 

concentration of the pollutant Co i.e.— — is proportional to 
V J at 

0 

— g— (C 2 -C 2 ) where v^ is the exchange velocity, (Reible et 
“e 

a7.,1983). In such a case the equation governing the concentration 
Cg of the pollutant can be written as. 
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+ u 


ac. 
3 aF 




- - ISjC, - r (C - C ) 
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C 2 (x,y,z,t) =0 at t = 0 


(4.59) 

(4.60) 


C 2 (x,y,z,t) = Cj(x,y,z,t) at x = x^ 


(4.61) 


CgCx.y.z.t) 
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as y 


± 00 


(4.62) 


for X <x„, and z i E . 

The concentration in region II may also change through 
inter-layer exchange characterized by an exchange velocity but we 
assume that there is no transfer of pollutants from the greenbelt, 
back to the adjacent thin layer as greenbelt is assumed to be a 
sink. Hence the equation governing the concentration in the 

greenbelt region remains the same as in the case 4.1. Further the 
equations for and also remain same as in the case 4.1. 

The equations ( 4 . 59 )- (4 . 62 ) can be written in dimensionless 
form as, (dropping bars for convenience) 
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C 2 (x,y,z,t) =0 at t = 0 


(4.64) 


C 2 (x,y,z,t) = Cj(x,y,z,t) at x = x^ 


(4.65) 



the solutions for and are given as follows whereas the 

solutions of Cj and C 2 are same as given above in section 4.1. 

4.3.1 CASE I WHEN FLUX IS INSTANTANEOUS AT THE SOURCE 
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4.3.4 DISCUSSION 

Thr pxj.rr.s.mr.ns ahtained above in this case have been 
compui<‘fi with t tic Kame set of parameter values as used for the 
case of 4.1. ft is found that in all the corresponding figures 
(not all figure.H shown) as to the previous section, the 
concentration behavior i .s similar (No significant difference in 
concentration has be<>n noted). In particular, in Figures (4.18) - 
(4.20) thf* f'oneent rat ion distribution of the pollutant at the 
receptor x = 0.1 and / = 0.02 for optimization of greenbelt 

location and the effeet of inter- layer exchange velocity v for 

© 

the insiantanenu.s flux at the source for three heights of the 
source i.e. , h > h , h = h^ and h^ < h. are shown. In this case 
also similar result.s have been found for the optimization of the 
greenbelt location nn in I igure.s (4.2) - (4.4). As the inter-layer 
exchange velocity inerease.s the concentration of the pollutant 
decreases at the receptor that means more pollutants from the 
region III ari' trapped in the greenbelt and since the transfer of 
the pollutants. from the region II, back to region III is 
neglected, the I'oncent rat i on reduces considerably. 

The (‘oncentrat it>n behavior for the case of constant and step 
function type flux is similar and hence not shown here. The 
effects of other removal parameters are also same as in the 
previous section. 

4.4 MAIN BESULTS 

In this chapter the effect of greenbelt on the reduction of 
concentration of pollutant has been studied when greenbelt is 
located between the source of the pollutant and the receptor. The 
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a’«' nblHincd by analyzing the model, 

( 1 ) The Ir.nK'itu.iinal rn„r r„t ra f i on of the pollutant along the wind 

direction always docrenses in presence of greenbelt and this 
decrease is further enhanced as the magnitude of depletion 
coeffJCKUit parameter due to greenbelt and its size (height and 

width) iiH'reasi 

(2) Similarly ver’nu.l roncentrat i on at a given location, after 

the greenbc'lt , det r« ases as the value of aforesaid parameters 

increase. 

(3) It has lieen .shown by using numerical optimization that the 
reduction in jn Mutant concent rat i on at the receptor is maximum 
under the followinK eonditions, 

(i) When the snuree height js greater than the greenbelt height, 

the gr(*enh('it should be loraled as close as possible to the 

receptor . 

(ii) When th«> sourci' height is less than or equal to the greenbelt 

height, the greenbelt should be located as close as possible to 

the source. 

(4) It has also been noted in both the above cases that for 
maximum ('oni'ent rat i on reduction at the receptor, the height of the 
receptor should b<' a.s small as possible than the greenbelt height. 

(5) The reduction in concentration at the receptor is larger when 
the deposition velocity at the top of the greenbelt or the 
inter-layer exchange velocity in the layer adjacent to the top 
of the greenbelt incr<iases. 
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FIG. 4.1 DISPERSION OF AIR POLLUTANT IN PRESENCE OF 
GREENBELT. 
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DIMENSIONLESS DISTANCE OF GREENBELT'FROM THE SOUiCE(xj 

FIG.4.2 CONCENTRATION DISTRIBUTION AT x=0.1 (RECEPTOR) FOR! 
hs>hg , WHEN FLUX IS INSTANTANEOUS. 


if 





DIMENSIONLESS DISTANCE OF GREENBELT FROM THi#1;e(x,) 
FIG.4.3 CONCENTRATION DISTRIBUTION AT x = ai (R^CHt‘''‘^) TOR 

h,=h,, WHEN FLUX IS instantaneous. 



z 

o 

k 

cr 


z 

lU 

(J 


8 

{/) 

U) 

LU 


Z 

g 

(/) 

z 

LU 


Q 



dimensionless distance of greenbelt from the source (x^) 

FIG.4.4 concentration DISTRIBUTION AT x=0.1 (RECEPTOR) FOR 
hs<hg,WHEN FLUX IS INSTANTANEOUS. 
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DIMENSIONLESS DISTANCE OF GREENBELT FROM THE SOURCE (x,) 

FIG.4.5 CONCENTRATION DISTRIBUTION AT x=0.1 (RECEPTOR) FOR 
hs>hg, WHEN FLUX IS CONSTANT. 




DIMENSIONLESS CONCENTRATION 


136 



DIMENSIONLESS DISTANCE OF GREENBELT FROM THE SOURCE (x,) 

FIG.4.6 CONCENTRATION DISTRIBUTION AT x = 0.1 (RECEPTOR) FOR 
hs = hg,WHEN FLUX IS CONSTANT. 
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DIMENSIONLESS DISTANCE OF GREENBELT FROM THE SOURCE (x,) 

FI6.4.7 CONCENTRATION DISTRIBUTION AT x=0.1 (RECEPTOR) FOR 
hs<hg,WHEN FLUX IS CONSTANT. 




DIMENSIONLESS CONCENTRATION 


0.35 
0.30 
0.25 

0.20 
0.15 
0.10 
0.05 
0 

DIMENSIONLESS DISTANCE OF GREENBELT FROM THE SOURCE (x^ 



FIG.4.8 CONCENTRATION DISTRIBUTION 
hs>hg,WHEN FLUX IS STEP 


AT x=0.1 (RECEPTOR) FOR 
FUNCTION TYPE. 
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DIMENSIONLESS DISTANCE OF GREENBELT FROM THE SOURCE(;( 

FIG.4.9 CONCENTRATION DISTRIBUTION AT x=0.1 ( RECEPTOR) 
FOR hs= hg , WHEN FLUX IS STEP FUNCTION TYPE . 
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DIMENSIONLESS DISTANCE OF 6REENBELT FROM THE S0L)RCE(x,) 

FIGA.IO CONCENTRATION DISTRIBUTION AT x = 0.1 (RECEPTOR) FOR 
hs<hg .WHEN FLUX IS STEP FUNCTION TYPE. 
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DIMENSrONLESS CONCENTRATION 

FIG.4.12 VERTICAL CONCENTRATION PROFILES AT x = 0.1 

(RECEPTOR) FOR hs = hg ,WHEN FLUX IS CONSTANT. 
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DIMENSIONLESS CONCENTRATION 

FIG. 4.13 VERTICAL CONCENTRATION PROFILES AT x=0.1 

(RECEPTOR) FOR hs<hg,WHEN FLUX B CONSTANT. 
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DIMENSIONLESS DOWNWIND DISTANCE 


FIG. 4.14 CONCENTRATION DECREASE IN GREENBELT FOB 
CONSTANT FLUX. 
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DIMENSIONLESS CONCENTRATION 

FIG,4.15 EFFECT OF X ON CONCENTRATION DISTRIBUTION 
AT xrO.1 (RECEPTOR) 
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DIMENSIONLESS CONCENTRATION 

FIG.4.16 EFFECT OF Nj.THE DEPOSITION VELOCITY AT THE TOP 
OF THE GREENBELT, ON CONCENTRATION DISTRIBUTION 
AT X =0.1 (RECEPTOR) 
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DIMENSIONLESS CONCENTRATION 

FIG.4.17 EFFECT OF INCREASING RECEPTOR LOCATION FR01 
THE GREENBELT AT x,= 0.04. 
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DIMENSIONLESS DISTANCE OF GREENBELT FROM THE SOURCE(x,) 


FIG.4.18 CONCENTRATION DISTRIBUTION AT x=0.1 (RECEPTOR) FOR 
hs>hg,WHEN FLUX IS INSTANTANEOUS. 
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DIMENSIONLESS DISTANCE OF GREENBELT FROM THE SOURCE(X|) 

FIG.4.19 CONCENTRATfON DISTRIBUTION AT x = 0.1 (RECEPTOR) F(^ 
hs = hg,WHEN FLUX IS INSTANTANEOUS. 
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DIMENSIONLESS DISTANCE OF GREENBELT FROM THE SOURCE (x,) 

FIG.4,20 CONCENTRATION DISTRIBUTION AT x=0.1 (RECEPTOR) FOR 
hs<hg .WHEN FLUX IS INSTANTANEOUS. 




CHAPTER V 


Effect of Precipitation Scavenging on Dispersion of Air Pollutant 

IN THE Atmosphere 

5.0 INTRODUCTION 

The study of various kinds of removal mechanisms which occur 
in the environment such as scavenging by rain droplets, wet and 
dry deposition of both particulate and gaseous pollutants is very 
important in view of the fact that our environment is continuously 
stressed by different kinds of air pollutants/toxicants emitted 
from different sources. One such important mechanism for removal 
of pollutants from the atmosphere is due to precipitation 
scavenging , a process useful in cleaning the environment of 
polluted cities in India during monsoon season. The study of 
scavenging of pollutant by precipitation has received considerable 
attention during the last few decades, e.g., Hales(1972 , 1989) , 
Hales et ai. (1973), SI inn( 1974 , 1977) . Peters( 1976) , Hill and 
AdamowiczC 1977 ) , AdamowiczC1979) , Overton et ai.(1979), 
Fisher ( 1982) , Scot t ( 1982a , b) , Levine and Schwartz(1982) , Shukla et 
ai.(1982), Chang(1984), Kumar ( 1985 , 1986 , 1989) , Shukla and 

Chauhan(1987) . Luecken et a/. (1989). In particular, Hales(1972) 
presented some fundamentals for the general analysis of 
atmospheric gas scavenging by rain and pointed out that the 
process of trace gas scavenging is reversible in nature and 
phenomenon of absorption and desorption may cause a redistribution 
of pollutants in the atmosphere. A simpler linearized model for 
nrftr? i 1 +wo uinohnnt nf ifaseous constituents from industrial 
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plumes by taking into account the reversibility of adsorption 
process has been presented by Hales g t a7#(l973) assuming the rain 
to be composed of spherical, non-interacting drops of static size 
distribution. An analytical study for the redistribution of gas 
plume caused by reversible washout has been investigated by 
Slinn(1974). The problem of irreversible scavenging for HNO„ vapor 
and r 'culation of washout coefficients for removal of HNO,. from 
the ^‘mosphere has been considered by Levine and Schwartz(1982) . 
Kumari.1985) studied the simultaneous process of trace gas removal 
from the atmosphere and absorption of these gases in rain drops 
using an Eulerian model but his model does not take into account 
the wind velocity and presence of a time dependent source. 

In the present investigation we, therefore, study the effect 
of precipitation scavenging on the unsteady state dispersion of a 
reactive air pollutant emitted from a time dependent point source, 
assuming the uniform distribution of rain droplets in the 
atmosphere which absorb the pollutant and remove it by their fall 
on the ground with a constant velocity. We consider that the 
concentration of the pollutant in the gaseous phase is governed 
by the unsteady state three dimensional convective diffusion 
equation with removal terms while the concentration of the absorb 
pollutant in the droplets is governed by a simplified equation 
which neglects diffusion term but takes into account the falling 
of rain droplets on the ground. These equations are then solved 
using appropriate initial and boundary conditions to study the 
effect of precipitation scavenging on the dispersion of pollutant 
in the atmosphere. 
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analysis are: (i) Instantaneous, and (ii) Constant. 


5.1 mathematical model 

Consider the unsteady state dispersion of a reactive air 
pollutant emitted from a time dependent point source of strength 
Q(t) located at height h^ and the washout takes place below the 
cloud base situated at height z = H. Here we consider two phases 
in the atmosphere viz., the pollutant in the gaseous phase and 
absorbed pollutant phase in the droplets. 

The differential equations governing the concentration 
C(x,y,z,t) of the pollutant in the gaseous phase is written as. 


ac ^ ac _ „ a^c 
al ^ ax ■ ^y rr 


K 


a^c 


ay 


^ az^ 


Ck + k^)C + pC^ 

g a 


(5.1) 


where x, y, z are cartesian co-ordinates, t is the time, x-axis is 
taken along the wind direction, z-axis vertically upwards, U is 
the mean wind velocity, K and K are diffusivities in y and z 

y 2 

directions respectively, k is the first order removal rate 
coefficient of the pollutant due to some natural mechanism and 
is the depletion rate coefficient due to precipitation scavenging 
by rain which is zero if their is no rain and p is the coefficient 
of rate of desorption of the absorbed pollutant C^ in the droplet 
phase . 

We prescribe the following initial and boundary conditions 
for equation (5.1), 

(5.2) 


C(x,y,z,t) = 0 


at t =0 


C(x,y,z.t) = 6(y) 6(z-h„) at x = 0 

u s 


(5.3) 


(5.4) 
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dC 

dz 

= 0 

at z = H 




(5 

.6) 


The 

f o 1 lowing forms 

of t 

ime 

dependent 

flux Q(t) 

at 

the 

point 

source are considered in 

the 

analysis , 



(i) 

Flux 

is Instantaneous, 

Q(t) 

= 

5(t) 

(5, 

.7) 

(ii) 

Flux 

is Constant, Q( 

t) = 



(5 

.8) 


Now 

we proceed to 

write 

the 

equation 

governing 

the 


concentration of the absorbed pollutant in the droplet phase. 
Since the concentration of the absorbed pollutant in the droplets 
is directly influenced by the concentration of the pollutant in 
the gaseous phase, it is necessary to consider simultaneously the 
coupled process of pollutant depletion in gaseous phase and 
absorption in the droplet ph se of the pollutant species, (Kumar, 
1985). 

The equation governing the concentration C (x.y.Zjt) of the 

£l 

absorbed pollutant in the droplet phase is written as, 


^ = > c - P 


(5.9) 


with initial and boundary conditions as follows, 

C (x,y,z,t) =0 at t = 0 (5.10) 

CL 

Cg^(x,y,z,t) =0 at z = H (5.11) 

where w is the terminal velocity of rain droplets falling on the 

ground, which depends on its radius r but in the following 

analysis we assume the terminal velocity as constant and rain 
droplets to be uniformly distributed. The absorption rate 
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removal rate coefficient, are also assumed to be constant. 

In equation (5.1) represents the source term due to 

desorption of the absorbed pollutant. It may be noted that ^ rj 
and u/k « 1. therefore, we neglect the desorption of the absorbed 

s 

pollutant (i.e., is assumed negligible). 

Introduce the following dimensionless quantities. 


t = 




X 


K X 
z 

^ 

U 


y 


~W~ ’ ^s 


s 

TT 


P _ U 

■ “57 


W = 


Hw 


C, = 


H Y 


U H" , _ Q(t) 

Ca- Q(t) - 


„2 

- _ H n 

n - -V 

Z 


Q -’o 


The system of equations (5.1-5.11) can be written in 
dimensionless form as, (dropping bars for convenience) 


dc ac _ a^c a^c _ . + „ )c 


C(x,y,z,t) = 0 


at t = 0 


(5.12) 

(5.13) 


C(x,y,z,t) = Q(t) 5(y) 5(z-h ) at x = 0 


C(x,y,z,t) =0 as y — > ± <s> 


dC 

dz 

dC 

dz 


NC 


= 0 


at z = 0 

at z = 1 


(5.14) 

(5.15) 

(5.16) 

(5.17) 


ac ac 

“ W = y C - T] C 


at 


az 


Cg^(x,y , z ,t)=0 at t=0 


(5.18) 

(5.19) 


Ca(x,y,z,t) 


0 at z 


1 


(5.20) 



156 


where a - — < « 
z 


' ^ " "IT > ^ = — 4~ 
z z 


(i) Flux is Instantaneous, Q(t) = 5(t) 

(ii) Flux is Constant, Q(t) = 1 


(5.21) 

(5.22) 


5,2 CONCENTRATION DISTRIBUTIONS 

The solutions for concentration distribution of pollutant, C 

and C are obtained by solving equations (5.12) - (5.20) and 
a 

written in each cases as follows, 


5.2.1 CASE I WHEN FLUX IS INSTANTANEOUS AT THE SOURCE 


C(x,y,z,t) = 


■y^/4px 


9 ^ 

2s_ > n s , 


exp^ - (a-t-ag+X^)xj -p — ^ 6 (t-x) 


-y /4px « . 2 'i ^n^^s^ 

C_(x,y,z,t) = Q^y ) exp ( - (a+a,,,+X„)x) f^(x,z,t) 

^ ° A4npx ^ ^ ^ n 

(5.24) 

where 


(5.23) 


F„(z) = cos X z + — r — sin X„z 
n n X n 

n 


= (cos X z + — r — sin X„z) dz 
n J n n 


X^® are the roots of the following transcendental equation 


X^ ^ ^ =1,2,3, 


f^(x,z,t) = - f^j (x , z , t) l^cosX^ + N — — J 

sin X^Za 

+ f„3(x,z,t)[oos J.^2 * N J 

+ f , z , t ) In cos x^z - x^zj 
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= e-"'^-''>cos(x„w(t-x - -L^)) H(t - X - 13 ]H(t - 13] 

= x^e-'’‘‘-'‘>sin(x^w(t-x - l 3 ,] H(t - x -l 3 ]H(t - l 3 ] 

fngCx.z.t) = e“’^^^"^^cos[x^w{t - x)] H(t-x) 

f ,(x,z,t) = -4- e"'^^^“^^sinrx„w(t - x)! H(t-x) 
n4 (. n J 

and 6(t-x) is the Dirac delta function defined as, (Carslaw and 

Jaeger, 1941) 

6(t-x) =0, t < X 

= 1/g, X i t 5 x+e 
= 0, t > x+c 

5.2.2 CASE II WHEN FLUX IS CONSTANT AT THE SOURCE 


-y /4px « , 2 X 

C(x,y,z,t) = ) expl - (<x+a^+X„)x] H(t“X) 

44npx ^ g n 7 1 ^^ n 

(5.25) 

_y2/4|5x ® 

C„(x,y,z.t) = y- > expf - <a+a^+X„)x] F„(z)g„(x , z . t ) 

J4npx 4:, ^ ^ n 

^ (5.26) 


where 


s inX. 


^ tJt X 

g^(x,z,t) = - g^j (x , z , t) l^cosx^ + N — ^ 

+ [ “ N cosXj^ + XjjSinXj^^ 


sin X_z^ 


^ Kill 

+ gj^3(x,z,t) l^cos Xj^z + N — J 

+ gj^4 (x,z , t) I^N cos X^z - X^sin X^zj 
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g , (x , z , t ) 
®nl 




nC i~z) 

-rjCt-iJJ- oos(argl) + sin(argl)j-tTi e ” 

i 2 2 

1 *1^ w 

, H(t - X- H(t - ‘^) 

_ r|( 1-z) 

1 sin(argl) “X^w cos(argl )j +x^w e 


2^2 2 
n +X w 
' n 

l-z> 


8„3(x,z.t) = [ 


X H(t- X- H(t 

-r)( t-x) ^ 


1-z 

w 


rj cos(arg2) + X^w sin(arg2)j +r} 


2~T2 2 

n *x„ w 


H(t-x) 


( X , Z I t ) 


= - f- 


-rjC l-x ) 


^ -r} sin(arg2) - X^w cos(arg2)j +X^W n 
S'! 7s 7i H( t-X ) 


^ 2 2 


argi = X„wCt- X - -i^). arg2 = X w(t - x) 

n w n 

and H(t-x) is the Heaviside function defined as.CCarslaw and 
Jaeger , 1941 ) . 


H(t-x) =0, t X 
t 


X + G 

= 1 , t Z: x + e 


X < t < X+G 


5.3 DISCUSSION AND RESULTS 

To study the effect of precipitation scavenging on the 
dispersion of air pollutant in the atmosphere, the expressions for 
concentration distribution of the pollutant in the gaseous phase 
and that of the absorbed pollutant in the droplets phase have be 
computed and displayed graphically in Figures 5.1 5.6. 
dimensionless parameter values used in the computation 

1.0, a = 2.00, = 20.0. p = 10.0, N = 0.04, w = 1000, h 

8 


0.2 
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When flux is instantaneous at the source, the central line 
concentrations of the; pollutant in gaseous phase versus downwind 
distance with and without rain are shown in Figure 5.1 for 
different values of t. It is noted that the concentration of the 
pollutant decreases as time increases in both the cases and the 
concentration of the pollutant decreases considerably due to 
precipitation scavenging by rain. 

When flux is constant at the source, the central line 
concentrations of the pollutant in gaseous phase versus downwind 
distance with and without rain are shown in Figure 5.2 for 
different values of t. It is found that the concentration of the 
pollutant increases with increase in time and reaches its steady 
state value as t — > ®. It is noted that the concentration of the 
pollutant reduces considerably in this case also. 

The concentration of the absorbed pollutant in rain 
droplets along the center line (0,0,h ) with downwind distance for 
both instantaneous and constant source is shown in Figures 5.3 and 
5.4 respectively for different values of t. It is observed that 
this concentration decreases as time increases for the case of 
instantaneous flux and increases for the case of constant flux and 
reaches its steady state value as t — > oo. 

In Figure 5.5 the concentrations of the pollutant in gaseous 
phase with and without rain are plotted with vertical height at a 
particular time t = 0.06 for three values of x = 0.03, 0.04, 0.05 
for the case when flux is constant at the source. The 
concentration of the pollutant decreases as x increases at a time 
t > X. It is also seen that the concentration of the pollutant 
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In Figure 5.65 the variation of concentration of the absorbed 
pollutant in rain droplet phase with vertical height is also 
depicted for different values of x at t=0.06. It is noted that 
this concentration increases with decrease in height but decreases 
as X increases for t > x. 

From the above analysis, it may therefore be concluded that 
precipitation scavenging helps in cleaning the atmosphere in those 
regions of the world where precipitation takes place very 
frequently . 
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DIMENSIONLESS DOWNWIND DISTANCE 

- IG. 5.1 CONCENTRATION -DISTANCE *thE^S(^cI 

PHASE FOR INSTANTANEOUS FLUX AT THE 50UK(-t. 
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DIMENSIONLESS DOWNWIND DISTANCE 

FIG. 5.2 CONCENTRATION -DISTANCE PRO'^“-ES IN 

PHASE FOR CONSTANT FLUX AT THE SOUK 



DIMENSIONLESS CONCENTRATION 



DIMENSIONLESS DOWNWIND DISTA 

FIG. 5.3 CONCENTRATION-DISTANCE PROFILES 
PHASE FOR INSTANTANEOUS FLUX A 
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DIMENSIONLESS DOWNWIND DISTANCE 

FIG. 5.4 CONCENTRATION-DISTANCE 

DROPLET PHASE FOR CONSTANT FLUX AT TMt 

SOURCE . 




DIMENSIONLESS VERTICAL HEIGHT 



DIMENSIONLESS CONCENTRATION 

;,5.5 VERTICAL CONCENTRATION 

PHASE FOR CONSTANT FLUX AT THE SOURCE, 
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CHAPTER VI 


Effect of Precipitation Scavenging on Dispersion of Air Pollutant 

IN a two Patch Environment 


6.0 INTRODUCTION 

Since the precipitation scavenging is one of the important 
mechanism to remove the pollutants from the atmosphere, it has 
received considerable attention during the last few years as 
described in chapter V, [Hales( 1972) , Scriven and Fisher(1975a,b) , 
Slinn( 1974 ), Kumar(I985), .Shukla et ai. (1991a)], where we have 
discussed the effect of precipitation scavenging on the unsteady 
state dispersion of a reactive air pollutant emitted from a time 
dependent point source without considering patchiness in the 
atmosphere which arises naturally in the atmosphere with different 
topographical and meteorological conditions at different places. 

In this chapter we, therefore, study the effect of 
precipitation scavenging on the unsteady state dispersion of a 
reactive air pollutant emitted from a time dependent point source 
in a two patch environment assuming that rain droplets, uniformly 
distributed in the second patch, absorb the pollutant and remove 
it by their fall on the ground. ¥e consider that the concentration 
of the pollutant in both the patches is governed by the unsteady 
state three dimensional convective diffusion equations with 
different removal and other parameters . However , the concent ivat ion 
of the absorbed pollutant in the droplets in the second patch is 


Part of this chapter is 
International Conference 
Atmosphere-^Surface Exchange 
Washington, U.S.A. 


to appear in the Proceedings of Fifth 
on Precipitation Scavenging and 
Processes, 15-19 July, 1991, Richland, 
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governed by a simplified equation which neglects diffusion term 
but takes into areount the falling of rain droplets on the ground. 
To determine the eoneent rat ion distributions, in the two patches 
i.e., without rain and with rain, these equat ions are sol ved us ing 
appropriate initial and boundary conditions, with a view to study 
the effect of precipitation scavenging in cleaning the 
env i ronmcnt . 

As before, we consider the time dependent flux at the source 
to be of the forms: (i) Instantaneous, andCii) Constant. 

6.1 MATHEMATICAL MODEL 

Consider the unsteady state dispersion of a reactive air 

pollutant emitted from a time dependent point source in a two 

patch environment. It is assumed that the point source of strength 

Q(t) is located in the first patch at a height h whereas the rain 

s 

takes place in the second patch below the cloud base located at 
height z = H. (see Figure 6.1) 

The differential equations governing the concentration 
distribution in both the patches are written as, 


Patch I (0 i X <Xj , 0 :S z s H) 

The concentration Cj(x,y, 2 ,t) of the pollutant in the gaseous 
phase in the first patch is governed by, 


a^c, 


a^c. 




ac^ ac 

aT” “i oE jt X X.X 

where x, y, z are cartesian coordinates, 
taken along the wind direction, z-axis 
the mean wind velocity, and are 

y and z directions respectively and 


( 6 , 1 ) 

t is the time, x-axis is 
vertically upwards, is 
constant dif fusivities in 
kj is the removal rate 
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coefficient of the rwllutant due to some natural mechanism. 

Thf‘ initial and boundary conditions for equation (6.1) are, 

C j ( X , y , 7. . t ) = 0 a t t = 0 (62) 

Cj(x,y,7. .t) = ^>(y> 5(z-hg) at X = 0 (6.3) 

Cj(x,y,z,t) = 0 as y — > ± od (6.4) 

ac 

^zl aF“ ■ '"’dl^l z = 0 (6.5) 


aCj 

^ =0 at z = H (6.6) 

where v^j is the deposition velocity of the pollutant in the first 
patch at the ground z = 0 and Si.) is the Dirac delta function. 

The following forms of time dependent flux are prescribed at 
the source as in the previous chapter, 

(i) Flux is Instantaneous, Q(t) = Q^5(t) 

(ii) Flux is Constant, Q(t) = Q 

c 

Patch II(x 2 :Xj, OsziH) 

In this patch, where rain takes place, we assume that 
uniformly distributed rain droplets absorb the pollutant and 
remove it by their fall on the ground. The concentration of the 
pollutant in the gaseous phase is governed by the unsteady state 
three dimensional convective diffusion equation with removal and 
other parameters. 

The concentration C 2 (x.y,z,t) of the same pollutant in the 
second patch is governed by, 



ac,, 


,ii ' . , 




^ ' "'2"kg>Cj 
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(6.7) 


where u., Ks tlH me-HH wind velocity, ^nd £,2 constant 

diffusjvit M'N in > and ? directions respectively and is the 
natural n^rrxival rate coefficient and k is the depletion rate 
coeffKienf of tin pollutant due to precipitation scavenging by 
ram. 

The initial and tmundary conditions for equation (6.7) are, 


Cg ( X , y , z , t ) = 0 a t t = 0 ( 6 . 8 ) 

C 2 ( X . y , z . t ) = f j { X , y . z . 1 ) at = x ^ (6.9) 

( X I y f z * t ) “ P as y “*.4 ( 6 , 10 ) 


dCy 

^z2 ‘ 3z ” '*'d2^'2 a t z s: 0 (6.11) 

ac„ 

= 0 at z = H (6.12) 

where ts the deposition velocity of the pollutant in the 

second patch at the ground z = 0. 

The concentration C (x,y,z,t) of the absorbed pollutant in 

a 

the droplet phase is assumed to be governed by the following 
equation , 


dc^ ac^ 

^ - w ^ = y C2 - 17 

with initial and boundary conditions, 


(6.13) 


C..(x,y.z,t) =0 at t = 0 

=0 at z = H 

where w, assumed to be constant, is the terminal 
droplets falling on the ground, y (s k^) is the 


(6.14) 

(6.15) 
velocity of rain 
rate coefficient 
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of incrcaKH of ( and r/ is its removal rate coefficient in the 
env i ronmrnt . 


Introduce the following dimensionless quantities, 


K, t 

T - -V 

K„ X 

zo 


_ y 


z 

h 

T- g 

i — " ■ ♦ A 

2 

Hu 

0 

t y 

- -Y- . 

z = - 

H ■ 

E = — sp 

s H 

u. = 5(t) 

_ Q(t) 

If 

o 

Q K 

0 zo 

r = 

u h2 

0 

Cji s i =1 , 2 


2 ’ 

«cH 

1 

Qc 

„2 
u H 

7r - 0 n 

w = 

H w 


h2 

H 7 


H^n 

^a 5^ ^a' 

k ’ 
zo 

1 = 

K ’ 
zo 

n = 

O 

IS) 

where K and 

zo 

IJ are 

0 

reference diffusion 

coefficient and wind 


velocity respectively. 

The equations (6.1)~(6.15) are written in dimensionless form 
as, (dropping bars for convenience) 


Patch I (.0 s X < X,, Oszil) 

1 


dC . 8C . 

+ U, -5— i 

dt 1 dx 


d^C, 


d^c, 


P 


1 TT' * ^1 _2 


ay 


az 


- «iS 


Cj(x,y,z,t> =0 at t = 0 


( 6 . 16 ) 

(6.17) 


Cj^(x,y,z,t) = 9112 .. 5(y) 5(z-hg) at x = 0 

Cj^(x,y,z.t) =0 as y — > ± oo 

ac, 

= ^ 1^1 z = 0 

acj 

aT- = 0 at z = 1 

and 

(i) Flux is Instantaneous, Q(t) = Q^6(t) 

(ii) Flux is Constant, Q(t) = 1 


( 6 . 18 ) 

( 6 . 19 ) 

( 6 . 20 ) 

( 6 . 21 ) 
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Patch II (x 

it M j 

* s 2 s 

1 ) 

ac2 

_ + U2 jjjT 

- = P; 

^ ay" " '2 

az" 

C 2 ( X , y , 2 , t ) 

= 0 

at t = 0 


C 2 ( X , y , z . t ) 

= S 

( .X . > . 2 . t ) 

at 

C 2 ( X , y , z , t ) 

= 0 

as y — » ± 

OD 

dC 

= N C 

32 ' 


at z = 0 


3C^ 

— £ = 0 

dz 


at 2 = 1 


and 




aT'-'^aF" 

= y C 

'2 - " 


C„ ( X , y , z , t ) 

a. 

= 0 

at t = 0 


C^(x,y,z,t) 

= 0 

at z = 1 


where 





(oc2+ag)C2 (6.22) 

(6.23) 

“ (6.24) 

(6.25) 

(6.26) 

(6.27) 

(6.28) 

(6.29) 

(6.30) 



6.2 CONCENTRATION DISTRIBUTIONS 

The system of equations (6 . 16)-(6 . 30) is solved for different 
forms of flux as described earlier and solutions are obtained in 
each cases as follows, 


6.2.1 CASE I WHEN FLUX IS INSTANTANEOUS AT THE SOURCE 
Cj(x,y,z,t) = A p(x,y, 2 ) «(t - 



Q(x ,y ,z) 



(6.31) 


C 2 (x ,y , z , t ) 


X 


x-x 


(6.32) 
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Q - 7 

C ( X , y , z , t ) = " 

a 


u 


M ( X , y . z , t ) 


(6.33) 


where 


P(x,y,z) 


2 

-y /(4 j. ) 

e 1 


471—— X ^ 


rTr 1 1 ^ n 


^ — y /4^(x) 

Q(x ,y , z ) = — — y exp 


■/ 4n0 


4n0(x) n = l 


- r__L_ 


2 


-j(x-xp 


Q 


n 


V „...r 1 


m= 1 


F (h ) 1 

— S ^ J F_(z)G_(z)dz 


m 


m n 


0 


M(x,y ,z) 


► 2 , . „ , ® 
■y /4^(x) 


f^(x,z,t) 


^ 4 n ^ ( X ) n=l 




OD 


in= 1 0 


^1 

F^( 2 ) = COS X Z + — r — sin \z 
n n X„ n 

n 


^2 

G (z) = cos u z + sin uz 

n 


2 

Pn = (cos x_z + — sin X_z) dz 


0 


■n~ X„ n" 

n 


N 


r * * * 2 2 

= J (cos Pj^z + - 77 - sin /i^z) dz 


0 


u ' n 

^n 


Here x and n are the roots of the following transcendental 


n ''n 

equations respectively, 
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Xn tanX^^ = Nj 


"" ^’2' n = 1.2,3,... 


f ( X , z . t- ) = 
n 


- f^j(x.z,t>[^coaM^+ ~ sin^J 

+ ^2 cos^^+ sin/jJ 

r ^9 X 

+ f j^^aCx ,z , t) j^cos ^^z + _ sin 


^^v...,w (.,2 Af^z - sin 


+ f „ . ( X , z , t ) N„ cos 




, X X'”X 

f^j(x,z,t) = exll cos(argl) HU - jji- + — -1 

Ji 2 


f„2(x,z.t) = sin(ar«I) H^t ^ 




n 

, Z , t ) 

= exll 

f C X , z , 1 1 

exll 

' n 


f 

f ^ 

argl = fj^Yi 

i - 

(-U 

arg2 = u wj 

f 

t - \ 

r ’^l 


K 

'■ “1 

exll = expj 

[-.I 

[t - 

P(x) = 

X . + 

p2 

“1 

1 

“2 



r ^1 

x-x . 

t - 

+ 

1 


d 

to 

r* 1 

r ^1 



' + 

1 ’ 


U 2 JJ 


+ _i + 

‘1 “2 






x-x. 


u. 


( X - X J ) 
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6.2.2 CASE 
C j ( X , y . 7 . t ) 

C 2 ( X , y , z , t ) 

( X , y , z , t ) 
where 

N { X , y , z , t ) = 


( X , z , t ) = 


X 


n WHEN FLUX IS CONSTANT AT THE SOURCE 


U , 


I' ^ X , y . z ) H 




(6.34) 


u , 


N { X , y , 2 , t ) 


•y /4^(x) 


CD 


J 

/ 4np(x) n=l 

7 




(6.35) 



(6.36) 

r“2-^“g*''2''n 1 

1 U 2 J 

1 (x - X ) 

i J 

j gj^(x,z,t) 


w 2 

V ..xof - 1 1 ''-'‘’s’ f‘ 

I. u, J 


j < X , 2 , t ) 


"^2 

cosu + — sin/j 


J 


gj^2^5c,2.t) N2 cos^^+ sinfij 

f ^2 ^ 

s jj" sin 

gjj4 (x , z , t ) (^N2 cos fj^z - fj^ sin /j^zj 


ex 11 {-r| cos(argl) +gj^w sin(argl)] + ri e 


-T}(l-z)/w- 


~~T~ 5 T 

rj + w 


X H 




exll [-rj sin(argl) cos(argl)] +fi^w e 


-r} ( 1 -z ) /w-i 


, 2 ^ 2 2. 
"n'” "n “ ’ 




H 
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? . t ) 


g . (x 
®n4 


z , t ) 


cxin p cos{arg2) sin{arg2)] +p‘ 

vf, r ’'i 

2 

n * 



+ 


cxlli-p sjn(arg2) cos(arg2)]+p w' 


r ’‘i 


H t- 



+ 



6.3 DISCUSSION AND RESULTS 

To see the effect of removal parameters, the expressions for 
concentration of pollutant in both the patches and that of the 

absorbed pollutant in the droplets have been computed and shown in 
Figures 6. 2-6. 9. The parameters values used in the computation are 

Qq “ 1.0, ttj = 2.00, - 2.0, ag= 20.0, ~ ^2 ~ 10.0, u^ ” ^2 ~ 

1.0, = ^2 ~ ^1 “ ^ ^2 “ ^ ~ ^s ^ ^ " 

18.0, p = 0 . 01 . Xj = 0 . 04 , y = 0.0, e = 0.005. 

When flux is instantaneous at the source, the central line 
concentrations of the pollutant in gaseous phase with and without 
rain with downwind distance are shown in Figure 6.2 for different 
values of t. It is observed that the concentrations with and 

without precipitation scavenging decrease as time increases but 

this decrease is considerably enhanced due to precipitation 
scavenging. 

The central line concentrations of the pollutant in gaseous 
phase versus downwind distance with and without rain for constant 
source are shown in Figure 6.3 for different values of t. It is 
noted that the concentrations of pollutant with and without ram 
increase as time increases and reach their respective steady state 
values as t — » oo. It is also concluded that the concentration of 
the pollutant reduced considerably due to precipitation in the 

second patch i . e . , x > = 0 . 04 , 



Th(> ronrcritrat ion of the absorbed pollutant C in rain 

£L 

droplets alorif? the center line (0,0, h^) with downwind distance 
for both instantaneous and constant source is shown in Figure 6.4 
and 6..) for different values of t. It is observed that this 
concentration decreases as time increases for instantaneous source 
and increases with time for constant source and reaches its steady 
state value as t — > ®. 

In Figure 6.6 the concentrations of the pollutant in gaseous 
phase with and without rain for constant source are plotted with 
vertical height at a particular time t = 0.10 for three values of 
X = 0.07, 0.08, 0.09. The concentration of the pollutant decreases 
as X increase's at a time t > x. It is also seen that the 
concentration of the pollutant decreases due to precipitation 
scavenging by rain. 

In Figure 6.7 the variation of concentration of the absorbed 
pollutant in rain droplet phase with vertical height is shown for 
different values of x at t = 0.10. It is noted that this 
concentration increases with decrease in height but decreases as x 
increases for t > x . 

When Xj = 0.001, i.e. when the second patch where rain takes 
place is closer to the source, the vertical concentration profiles 
of the pollutant in gaseous phase with and without rain for 
constant source are shown in Figure 6.8 for different values of x 
= 0.03. 0.04, 0.05. It is found that this figure is approximately 
same to Figures 5.5 of the preceding chapter without patchiness in 
the atmosphere showing the importance of the location of the 
second patch from the source where precipitation is taking place. 
Similar result is also valid for absorbed pollutant in the droplet 
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ptmsr for rcnsUnl „our.-r, (See Figures 6.9 and 5.6). 

Finally, f ri.m the above analysis anri 

nai>sis and the results in chapter 

V. I I maj' he rn, .eluded that precipitation scavenging is a very 
i.portant ..eehan , .sn, for the removal of the pollutant from the 
atmosphere „ g.ven receptor point located in the second patch. 
However, the rate of reduction depends upon the location of the 
patch and prevailing meteorological conditions. It may also be 
noted that the reduction in pollutant concentration is larger at 
the receptor when precipitation patch is closer to the source. 
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FIG.6.1 DISPERSION OF AIR POLLUTANT AND ITS REMOVAL 
BY PRECIPITATION SCAVENGING IN THE SECOND 
PATCH . 


DIMENSIONLESS CONCENTRATION 
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DIMENSDNLESS DOWNWIND DISTANCE 

FIG. 6.2 CONCENTRATION-DISTANCE PROFILES IN GASEOUS 
PHASE FOR INSTANTANEOUS FLUX AT THE SOURCE. 




DIMENSIONLESS CONCENTRATION 


lej 



DIMENSIONLESS DOWNWIND DISTANCE 

FIG.6.3 CONCENTRATION-DISTANCE PROFILES IN GASEOUS 
PHASE FOR CONSTANT FLUX AT THE SOURCE. 






DIMENSIONLESS CONCENTRATION 
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DIMENSIONLESS DOWNWIND DISTANCE 

FIG.6.4 CONCENTRATION-DISTANCE PROFILES IN RAIN DROPLET 
PHASE FOR INSTANTANEOUS FLUX AT THE SOURCE. 




DIMENSIONLESS CONCENTRATION 
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DIMENSIONLESS DOWNWIND DISTANCE 

FIG.6.5 CONCENTRATION-DISTANCE PROFILES IN RAIN 
DROPLET PHASE FOR CONSTANT FLUX AT THE 
SOURCE. 




DIMENSIONLESS VERTICAL HEIGHT 


IS-I 



DIMENSIONLESS CONCENTRATION 


FIG. 6. 6 VERTICAL CONCENTRATION PROFILES IN GASEOUS 
PHASE FOR CONSTANT FLUX AT THE SOURCE. 





DIMENSIONLESS VERTICAL HEIGHT 
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FIG.6.7 VERTICAL CONCENTRATION PROFILES IN RAIN DROPLET 
PHASE FOR CONSTANT FLUX AT THE SOURCE. 




DIMENSIONLESS VERTICAL HEIGHT 
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DIMENSIONLESS CONCENTRATION 

FIG.6.8 VERTICAL CONCENTRATION PROFILES IN GASEOUS 
PHASE FOR CONSTANT FLUX AT THE SOURCE. 




DIMENSIONLESS VERTICAL HEIGHT 
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DIMENSIONLESS CONCENTRATION 

FIG,6.9 VERTICAL CONCENTRATION PROFILES IN RAIN 
DROPLET PHASE FOR CONSTANT FLUX AT THE 
SOURCE. 




CHAPTER VII 


An Ecological type Nonlinear Precipitation Model for Removal 
Mechanism of Air Pollutant Forming Secondary Species in the 

Atmosphere 


7.0 INTRODUCTION 

Irs the last tno chapters we ha\'e proposed linear mathematical 
models using unsteady state three dimensional convective diffusion 
equation to study the effect of precipitation scavenging due to 
ram, on the dispersion of air pollutant from a time dependent 
point source. Linear models related to precipitation scavenging 
due to rain, fog or snow and similar situations have also been 
studied by HaIes(I972), Slinn(1974), Fishert 1982) , Shukla et 
ai. (1982), Kumar(l985), Pandis and Seinfeld(1990) , etc. However, 
it may be pointed out that in real situations during 
precipi tat ion , the number density of rain, fog or snow droplets 
changes as the intensity of precipitation increases and this 
change in number density affects the interaction process and 
should, therefore, be taken into account in the model. The removal 
term, in general, due to precipitation is not only proportional to 
the concentration of the absorbed pollutant but is also 
proportional to the number density of rain droplets in the 
atmosphere making the phenomenon nonlinear. 

Keeping in view the above, in this chapter a nonlinear 
ecological type precipitation model for removal mechanism of both 
primary and secondary pollutants is proposed by taking into 
account the interaction of droplet phase with the pollutant phase. 
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( iv) combined phase of both the pollutants absorbed in droplets. 

Let C^(t) be the number density of droplets of radius r in 
the droplet phase, C(t) be the concentration of the primary 
pollutant, Cp(t) be the concentration of the secondary pollutant 
and ( t ) be the concentration of the absorbed pollutant in the 
droplet phase. 

We now proceed to formulate the nonlinear model governing the 
evolution of these phases using the concepts of ecology related to 
growth rate, death rate and interaction process, Maynard 
Smith(1974), Korraoncly ( 1 986) . The ideas presented here are similar 
to as used by Freedman and Shukla (1991) in another context. [See 
also Shukla et ai. C1991c)]. 

In the monsoon season it is generally observed that, during 

precipitation, rain starts with less intensity and later it 

becomes more intense and continues so for days. In such a case the 

number density of droplets increases as time increases. Similar 

situation arises when snow falls on the ground. Also in the case 

of foggy environment the number density of fog droplets increases 

from some smaller value to a higher value when fog becomes more 

dense. Further, the growth rate of droplet phase may also change 

because of deposition of droplets on the ground or on vegetative 

canopies due to gravity, and this rate may be proportional to the 

number density of droplets as well as to the concentration of 

pollutants in the droplet phase (i.e., 

Under these conditions and assuming the continuity of the 

process, the differential equations governing the number density 

C (t) may be written, using ecological , type ideas, as follows, 
r 
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(Maynard Smith, 1974) 


'*‘^r 

= q - r(C^)C^ 


(7.1) 


whoro nr^) := + „C^. and thus r(0) = r . 0, ric ) = n > 0 

^ fli 

and q IS the rate of increase of droplets in the atmosphere 
assumed to be constant. It may be pointed out here that in real 
situations this rate is time dependent and may be periodic but in 
this analysis q may be considered as the mean of this periodic 
growth rate. 

To write the equations for C(t), C (t). and C (t) it is 

P £L 

assumed that the pollutant is emitted in the atmosphere at a rate 

Q( t ) , with natural depletion rate and converted to secondary 

pollutant with a rate XC. Thus, the rate of increase of secondary 

pollutant is XC which may also get depleted with a rate 6,C . The 

1 P 

main mechanism, however, for removing pollutants (both primary and 
secondary) from the atmosphere is through the droplets falling on 
the ground and this removal rate is considered to be proportional 
to the number density of droplets as well as the concentration 
in the droplet phase (i.e., p C C ). Further, it is assumed that 

fit F fit 

both primary and secondary pollutants are absorbed by the droplet 
phase in direct proportion to the concentration of the respective 
pollutants and the number density of droplets (i.e., aCC^ and 
respectively). We also assume that the absorbed pollutant 
with concentration C in the droplet phase may be removed from the 

fit 

total environment directly with rate and some of it may 

reenter (i.e., 05 C and ) the atmosphere to increase the 

a a 1 a a 

growth rate of primary and secondary pollutants if the process is 


reversibl e , 
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With these considerations the differential equations for 
C(t), Cp(t) and C^(t) are written as, 


dC 

dt 


(6^ + X)C - aCC^ + + Q(t) 


dC 


dt 

dC 

£ 

WT 


P _ 


5.C + XC - a,C + 6,6 C 

iP Ipr laa 


+ a, C C - p C C 
aa r Ipr ara 


(7.2) 

(7.3) 

(7.4) 


for C (0) ^ 0, C(0) ^ 0, C (0) ^ k C^(0), C (0) ^ 0, where 

* a, r p 

is the natural depletion rate coefficient of the primary 
pollutant in the atmosphere, 

is the natural depletion rate coefficient of the secondary 
pollutant in the atmosphere, 

is the natural depletion rate coefficient of the absorbed 
pollutant in the atmosphere, 

X is the conversion rate coefficient of primary pollutant to 

secondary pollutant in the atmosphere, 

a. are the reaction rate coefficients of primary and secondary 

pollutants in droplets respectively, 

u ii. u ) is the removal rate coefficient of droplets from the 
a 

atmosphere, 

p is the removal rate coefficient of the absorbed pollutant 

a 

in droplets from the atmosphere, 

k ^ 0 is a constant determining the initial pollutant 

concentration in droplet phase at t = 0, (which is 
naturally dependent on C^(0) and note that for C^(0) = 0 , 
k has to be zero), 

0 i ^ 1 are the reversible rate coefficients. 
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We now analyze the model (7.1), (7.2). (7.3), (7.4) for Q(t) 

B 0 and Q(t) . Q (constant) by using stability theory, (Lasalle and 
Lefschetz, 1961). The case Q(t) e 0 corresponds to the situation 
when the pollutant is emitted in the atmosphere by an 

instantaneous source. However, the case Q(t) e Q represents the 
emission of pollutant in the atmosphere with a constant rate, for 
example, by stacks emitting continuously. 

7.2 ANALYSIS OF THE MODEL 

The analysis of the above mode] will be presented in two 
parts. We will first analyze the model assuming that there is no 
secondary pollutant in the atmosphere and later with the formation 
of secondary pollutant as well. 

7.2a ANALYSIS WITHOUT SECONDARY POLLUTANT 

In this case, in the absence of secondary pollutant, the 
equations (7.2), (7.3), and (7.4) reduce to the following form, 

^ = - 5 C - aCC^ + e5„C^ + Q(t) (7.5) 

dt 0 r a a 

dC 

-3-r^ = - 6 C + aCC - u C C (7.6) 

dt a a r a r a 

C(0) it 0, C (0) ^ k C (0) 
a r 

We now analyze the model (7.1), (7.5), (7.6) for Q(t) e 0 and 
Q(t) E Q (constant). 

7.2.1 CASE I Q(t) E 0 (Instantaneous Source ) 

In this case our model has only one non-negative equilibria 

in the C -C-C space .namely, E (— ^,0,0). To determine the 
r a 0 r^ 

local stability of the equilibria, we compute the corresponding 
variational matrix as follows. 
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Using Routh Hurwitz theorem, from the above it is clear that 

locally asymptotically stable in C -C-C plane. 
0 r a 

However, we can say much more about in the form of the 

following theorem, 

THEOREM 7.2.1 If C^(0) > 0, then E^ is globally asymptotically 
stable with respect to the non-negative octant. 

PROOF : First we note 

dC 

= Cr = Q - r(c^)c^ 

^ “ - ■'o^r 

Hence 1 im C (t) s — 5— 

. r r 

t > 00 0 

Now consider 

C + = -f 6 C + 6„C^(1 - 6) + u C C 1 

aLoaa araj 

i -fs^C + 6„C (1 - 6)1 

VO a a J 

i -5„(C + C ) 
m a 

where 6 = minr6 ,5(1- 0)1 

m V 0 a J 

-6 t 

Therefore, C(t) + C^(t) i [c(0) + e ® , and hence the 

system is dissipative. From the above it follows that 

lim C(t) = lim C (t) = 0. Hence in the limit C^(t) is given by 

t KX> t KX> 

the solution of C^(t) = q - r^C^.. Since 0^,(0) > 0, the theorem 
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foil ows . 

This theorem implies that if the pollutant is emitted in the 
atmosphere by an instantaneous source (C^(0) > 0), all the 

pollutant would be removed from the atmosphere by the 
precipitation scavenging with the number density of droplets 
remaining at its equilibrium. 

7.2.2 CASE II Q(t) s Q > 0 (Constant Source ) 

In this case our model has one non-negative equilibria 
E ( C^ , C , C^ ) , 

To prove that E exist, we show that the system of algebraic 
equations , 



function of C„ . 

a 



Further from equation (7.12), C — > 0 as C — ♦ 0, and C is also an 
increasing function of . Hence two isoclines (7.11) and (7.12) 
would intersect as shown in Figure 7.1, and give positive 
equilibrium, provided the line segment P.,P„ < P.P.,. where 

O 1 *3 

coordinate of P„ is (Q/6 , 0). This condition can be written from 

o m 

equation (7.11) and (7.12) as, 


(6 + 66 ) ( r 6 + l>Q) 
m a 0 m 

{6 ( r 5 + aq) + 6 j>Q) 

moo 0 


{6„(r 6^+ l>Q) + qu„5 } 
a o ra am 


aq5, 


m 


where 6 = min{6 ,6 (1 - 6)1, see lemma 7.2.3. 

m 0 a 

From Figure 7.1 it is noted that as u increases the equilibrium 

8 

values C and C both decrease and even may tend to zero if u is 
a ^ 

very large. 

Now the variational matrix corresponding to E is given by M 
as foil ows , 


M = 


- r(C„) 

a 


- aC 


aC - P.C 
8 8 


-(5 + aC ) 

0 r 


«C. 


-.Cr 

65 . 


-(5^ + 

a a r 


For E we now find a sufficient condition for this equilibrium to 
be locally asymptotically stable in the form of the following 
theorem , 


THEOREM 7.2.2 Let the following inequalities hold, 

2aC < r^+ 

0 8 8 

(i; - p )C < 5„(1 - 6) 
a r 8 


(7.13) 


then E is locally asymptotically stable. 
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PROOF : If inequalities (7.13) hold, then by Gershgorin’s theorem 
(Lancaster and Tismenetsky, 1985) all the eigen values of M have 
negative real parts, and the theorem follows. 

We now show that under certain conditions E is globally 
asymptotically stable. We first prove a lemma. 


Lemma 7.2.3. The region 


0 , 


(C^,C,C):0^C ^_S_;0^C + C i-S- 


m) 


where 5^ = min^S^, 5^(1 ~ 0)j , is a region of attraction, 
Proof : As before 


lim C (t) i -ii- 
r r 

00 0 


and C + C :S-5(C + C)+Q 

Si Wi sl 


Hence lim (C + C ) 

• £l 

t > 00 


Q 


, proving the lemma. 


m 


THEOREM 7.2.4 Let the following inequalities hold, 

" ‘'a’-f " ‘ '■'o " " ‘'A’ 


m 

mo a 0 r 


(05^ -H aC^) < (5^+ u^C^)(5^+«C^) 


(7.14a) 

(7.14b) 

(7.14c) 


then E is globally asymptotically stable. 

PROOF : Consider the following positive definite function about 

i. 


= i cc_- c_)2 + 4 (c - c)2 + i (c^- c^) 


V(C^, C, C^) - J ^ 


(7.15) 


The derivative of V with respect to t along solutions can be 



198 


written as 


V = (C^- C^)(q - rcc^)cj + (C - C)[ - 6^C - <,CC_. t 66 c + 


+ (C - C ) r - 5 + aCC - i; C C ] 

a at aa r araj 


(7.16) 


After some algebraic manipulations, it can be written as 


V = - (r^+ + ((aC - 


1(C - £ )(C - C ) 

a a r) r r a a 


- (5^+ aC^)(C - C) - aC(C - C)(C^- C^) 

■ 'V ^a>^ * '“a =a> 

which can further be written as the sum of three quadratic forms, 

^ = - 5 ^a^ ’ I ^2'V ^a^' 

- 2 “22'Ca- ^a’^ * “23"= - C>fCa- Ca> " I 


- i a„.(C - C)^ + a,.,(C - C)(C^- 5 ) - i- a, . (C - 8)^ (7.18) 

Zoo Id rr Zilrr 


where 

^11 

= r + 
0 


^12 

(aC ■ 

■ ‘'a'^a> ■ ‘''^r 



^22 

= (V 

"aCr' ■ 

®^23 

66 + 
a 




^33 

= (5 + 
0 

aC^) , 

^13 

- aC 



Thus a 

suf f icient 

condition for 

V to 

be negative definite 

is that 

2 

^12 

ail 

^22 ^ 

0 




(7.19a) 

2 

^13 

^11 

^33 ^ 

0 




(7.19b) 

2 

^23 

^22 

®^33 ^ 

0 




(7.19c) 


hold. However (7.14a) implies (7.19a), (7.14b) implies (7.19b) and 
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(7.14c) implies (7.19c). Hence V is a negative definite and so V 
is a Liapunov function with respect to E, whose domain contains 
Op proving the theorem. 

This theorem shows that provided the inequalities (7.14) 

hold, the system settles down to a steady state and as the removal 

rate coefficient increases the concentration of the pollutant 

in the atmosphere decreases and it may even tend to zero if i> is 

a 

very large. 

7.3a ANALYSIS WITH SECONDARY POLLUTANT 

Now we will analyze the model ( 7.1), (7.2), (7.3), and (7.4) 
considering the formation of secondary pollutant in the 
atmosphere as well. 

7.3.1 CASE I Q(t) H 0 (Instantaneous Source ) 


In this case our model (7.1 - 7.4) has only one non-negative 

equilibria in the C - C - C - C space, namely, E. (■ ■ -- ,0,0,0). The 

r p a 1 r 

variational matrix corresponding to E^ is given by as follows, 



To determine the local and global stability behavior, we propose 
the following theorems. 

THEOREM 7.3.1 If the following inequalities hold, 

(9 + 9j)5^ < 5^+ 


(7.20) 



200 

then Ej is locally asymptotically stable. 

PROOF : The proof follows similarly as in Theorem 7.2.2. 

THEOREM 7.3.2 If C^(0) > 0, then E^ is globally asymptotically 
stable with respect to the non“negative octant. 

PROOF : First we note, as before 
dC 

-dt^ = = q - r(C^)C^ 

^ ^ 

Hence lira C (t) i ^ 

+ r r 

t > CO 0 

Now consider 


C + C + C 
P a 


SC + 5.C + (1 - e-e,)+pcc 

0 Ip aa 1 ara 


- [ «„C * 6jC 


p " W- - s-ep 




where 5^ = minfs , 5,, S (1 - 0 - 0. )! 

n V 0 1 a 1 J 

r > 

Thus C(t) + C (t) + C (t) £ fcCO) + C (0) + C (0)1 e . 

p a 1, p a j 

hence the system is dissipative. From the above it follows that 

1 im C(t) = 1 im C (t) = 1 im C (t) = 0, Hence in the limit C Ct) 

t X» t KD ^ t KX) ^ ^ 


is given by the solution of C (t) = q - r C , since C (0) > 0, 
and the theorem follows. 

This theorem implies that in the case of instantaneous 
introduction of primary pollutant in the atmosphere and 
consequential formation of secondary pollutant, the process of 
precipitation scavenging will remove both primary and 
secondary pollutants from the atmosphere completely. 
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7.3.2 CASE II Q(t)sQ>o (Constant Source ) 

In this case our model again has only one non-negative 
equilibria E*(C*, C? C* . C* ) in the C^-C-Cp-C^ space. 

To prove that E exist, we show that the system of algebraic 
equat ions , 


q - r(C*)C* = 0 

a r 


(7.21) 


-(6^+ X)C - aC*C* + + Q = 0 


(7.22) 


6 C* + XC*- a.cV + 0 6 C* = 0 
ip Ipr laa 


(7.23) 


- 5 C* + aC*C* + a,C*C* - p C*C* = 0 
aa r Ipr ara 


(7.24) 


has a positive solution, 


Equation (7.21) gives, C = , 

^ r(C ) 
a 

substituting C^ in equation (7.22) and (7.24), we get 


(7.25) 


* Q + 06 C 
C* = ^ a 

q 

X + 6 + a j — 

° r(C^ 

a 


(7.26) 


^ C* f5„r(C*) + qy ' 
C = ^ I a a ^ aj 

aq 


J- C* 
a p 


(7.27) 


where Cp, is given by, after using C^, in equation (7.23) 

C = 

P 5 + « ^ 

^ r(C ) 

a 


(7.28) 


Q . 05^C^ 


X + 6 + a- 


r(C ) 
a 


and 
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From equation (7.28) it is noted that C* increases as C* 

P a 

increases. And 
Q X 




X + 5 + oc*— ^ 

0 r 


j|c jlj 

= C . (say) > 0, as C 

a 


We also note from equation (7.28) that when X, , and e are 
zero, Cp is zero. Further we note from equation (7.27) that C* 
varies with through the first term only and thus the 

equilibrium point, if it exists, will have similar behavior with 
respect to as in the case of (7.2.2) without secondary 

pol lutant . 

Now we will show that the two isoclines (7.26) & (7.27) 

9)C 3fc 

intersect in the interior of C - C plane. 

a 

From equation (7.26) it is seen that C increases as C increases 

ft 


and C 


Q 


X+ 5 + a— ^ 
0 r 


= C. (say) > 0 as C 

i- Sl 


0 


From equation (7.27) it is seen that as C* — > 0, C* — > - — ^ C*. 

a ’ a pi 

SfC ^ 

since C — > C , as C — » 0. We also note from equation (7.27) 
p p 1 a 

♦ 

that to lie C in +ve quadrant, the condition. 


6 r(C ) + qu 
a a 


sl > 

— •- > a 


jCp , must hold. 


Hence two isoclines would intersect as shown in Figure 7.2 and 

give positive equilibrium provided, the line segment R 2^3 ^ ^1^3 

where the coordinate of R., is (Q/5 , 0) and this condition can be 

o n 

written by the equations (7.26) & (7.27) as, 


+ i>Q> * 


(7.29) 
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where 


( 6 + 0 6 )(r 6 +pQ) 

n a on ^ 


^n ~ niin{5^, 6^, 5^(1 ~ 6 ~ 6j)}, see lemma 7.3.4. 

From Figure 7.2 it is noted that as t> increases the equilibrium 

Si 


level of C and 

C " 

a 

decreases as 

in Figure 

7 . 1 and 

thus 

the 

equ i ] i br i um 1 eve 1 

of 

* 

C decreases 

P 

since C 

a 

decreases 

(see 

Eq. 

7.28) and hence 

both 

primary and 

secondary 

po 1 1 utants 

will 

be 


removed and may even tend to zero if p is very large. 

Si 

Now the variational matrix corresponding to E is given by M 
as foil ows , 


M 


- r(C ) 

£L 


- aC 


— C 5 + X + aC ) 

0 r 


0 

0 


r* 

-a . L 
1 P 

* * * 
aC + a.C - P„C 
Ip a a 


aC 


-(5.+a,C ) 

1 1 r 

“l^r 


-.Cp 

05. 


^l^a 


-(5, + P,C ) 
a a r 


THEOREM 7.3.3 Let the following inequalities hold 
2(„C% ajC*) < r( C*) .• .,^c’ 

(p - < 8^(1 - e - ep 

Then E is locally asymptotically stable. 


(7.30) 


PROOF : The proof follows similarly as in Theorem 7.2.2. 

We now show that under certain conditions E is globally 
asymptotically stable. We first prove a lemma. 
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Lemma 7.3.4. The region 


^2 » |(C^ ,C ,Cp, C^) : 0 i i ; 0 i C + C^+ i -. | j 


where 6 = min 

n 


[^ 0 ' ^ 1 ’ - 6 - )j . is a region of attraction. 


Proof : As before 


1 im C ( t ) £ — ^ 

CO 0 


and C + C+ C ^-S(C+C+C)+Q 
pa n p a ^ 


Hence 1 im (C + C + C ) 

P a 


Q 


cx> 


proving the lemma. 


THEOREM 7.3.5 Let the following inequalities hold 


[Ca + Oj + + »C*)( 5 ^ + p^C*) (7.31a) 


n 


„V ^ _4_ ^ )C5^ . X . aC^) 


(7.31b) 


(05^ + aC* )2 < 4_(5^. X ^ aC* ) 


(7. 31c) 


( 6 , 6 ^ * CjC* )2 < 4-(6^* OjC* ) 


(7.31d) 


^ -4 < " «i<) 


X^ < — 5 — ( 6 . + a.C )($ + X + aC ) 

a 1 1 r 0 r 


(7.31e) 

(7. 31f ) 


then E is globally asymptotically stable. 


PROOF : Consider the following positive definite function about E 
VCC^.C. Cp.C^) = i (C^- c* )2 + 4 (C - c*)2, 1 (Cp- c*)2 


+ — CC - C*)^ 
2 '•a a^ 


(7.32) 

The derivative of V with respect to t along solutions can be 


written as 
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V = (C^- Cj,)(^q - r(C^)cJ + (C - C*)[ - (5^+ X)C - aCC^ 

+ 65 C + q] + (C - C*) r - 6 C + XC - a, C C + 0,5 C 1 
J PPV Ip Ipr. 1 a aj 

+ (C - C*)r - 5^C + aCC^+ a,CC - i^Ccl 

a ax aa rlpr aray 


(7.33) 


After some algebraic manipulations, it can be written as, 


V 


<)(Cr- cy * ((aC . o,jCp-„^C^) -pcj (C^- cjxc^- C* ) 


- (5 + X + aC*)(C - C*)^ - «C(C - C*)(C - C*) 
or r r 


- (5,+ a.C*)CC - C*)^ + (05 + aC*)(C - C*)(C - C*) 

llrpp a r aa 

■ P„C*)(C - C*)^ + (0,5^ + a,C*)(C - C*)(C- C*) 

aaraa la Iraapp 

+ X(C- C*)(C - C*) - a,C (C^- C*)(c - C*) - 
pp Iprrpp 

which can further be written as the sum of the quadratics as 

- 2 <>«a- - 5 


1 Nl'V C*)" - aj3(C^- C’)(C - C*) - 4 a33(C- C*)^ 


^ - ,(C.- c!)^ + a^.,(C„ - c!)(C - C*) - i a 33 (C - C*)^ 


^ “22" a "^a- 


‘23"'^a a' 




i “n<V " ^ 4 'V <><v i ^ 44 ‘S' 4^^ 


^ a^^CC - C*}^ + a3^(C - C*)(C^ " ^p^ " I ^44^S' 


34 


P P 
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where 

1 

2 

3 


^2 

= (aC + a. C 

1 P 

- P C ) - pC 
a a r 


^22 

2 

3 


^3 

= - «C 



^33 

2 

" J 

($^+ X + aC*) , 

^23 

= 6S + aC 
a r 



^44 

2 

(6j + a^C*) , 

^24 

= + a^C* 



^14 

= - 

“iS’ 

^34 

= X 


Thus a 

sufficient condition 

for V to 

be negative 

definite is that 

2 

^12 

af j 

^22 

< 0 



(7.34a) 

2 

^13 - 

aj j 

^33 

< 0 



(7.34b) 

2 

^23 

^22 

^33 

< 0 



(7.34c) 

2 

^24 

^22 

^44 

< 0 



(7.34d) 

2 

^14 - 

^11 

a . ^ 
44 

< 0 



(7.34e) 

2 

^34 " 

^33 

44 

< 0 



(7.34f ) 

hold. 

However 

(7.31a - 

7.31f ) 

implies (7 

.34a - 7.34f) 


respectively. Hence V is a negative definite and so V is a 
Liapunov function with respect to E , containing the domain Q 2 > 
proving the theorem. 

This theorem shows that if there is constant emission of the 
pollutant then even in that case both the pollutants can be 
removed from the atmosphere and this rate of removal would depend 
upon the rate of introduction of the pollutant and other 


parameters . 
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7.4 CONCLUSION 

We have proposed a nonlinear ecological type model for 
removal mechanism of both primary and secondary air pollutants in 
the atmosphere by the precipitation scavenging. Using stability 
analysis it has been shown that when the pollutant is emitted in 
the atmosphere by an instantaneous source (single introduction of 
the pollutant in the atmosphere), the pollutants from the 
atmosphere would be completely removed by the precipitation 
scavenging. When the pollutant is emitted at a constant rate in 
the atmosphere it is found that the pollutants may be removed from 
the atmosphere completely but it would depend upon the rate of 
emission and removal by precipitation and other processes. 

The proposed model and analysis presented here also suggests 
that to remove the pollutants from the atmosphere, external 
species may be introduced in the atmosphere which can interact 
with the pollutant and remove it by some suitable removal 
processes based upon the physical and chemical properties of the 
pollutants as well as that of externally introduced species. 
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